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CHAP2ER I 
DELCAIPTION OF OOLID HELI 


imre Ces Colids 

The plyyoicel properties of the solid phases of the “rare pas” 
elemento ave of particular interest because their etemic electrona occupy 
Suily closed shells, and the only chemicei binding foree ia the weak ven 
der Weals interaction. The piysical properties of these solide are thus 
in many reapects more open to anmlysie than most other eclids. Amather, 
the weeiness of the interparticle force givea rice to deMinite quantum 
mechenical effects in the motion of the atoms im the condenaed state. 

Dovke azul Jonee™ heve written an extensive review of the experi- 
menind and theoretical work on there golids, with emphasis on s0iid argon. 
they show the importance of the quantimemechenicel sero-polnt energy by 
disting the ratio of tero-point energy to usheaive energy fox the rare 
geees (listeé in order of decreasing etonic weight): xenon, 6.031; 
krypton, 0.0543 argon, 0.0973 neon, 0.321; helium four, 3.92. Nenee, 
ostinary clescical inttice dynemics applies to mencn, Krypton, and argoa, 
with omail quantum-emechanical corrections. Neon can be treated alequntely 
by including anharmonic terme in the interatomic potential. However, 


helium is of such low atomic mages that more power techniques must 





i & Dobbs and C. ©. Jones, Theory and Properties of Colid 
Argon. Geports on Progress in Physics, Vol. EX, 526 (1957). 
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De empioyed.” Pucther, there i9 alse an isotope of ataxic weight three, 
vbich has spin 1/2 ang obeys the Pomal-Direc statistics. 

this theaia ig Givected to the task of deeevibing the growd state 
oe He? in @ way vhich includes quantum effects in the basic formelian. 
ane deseription of the ground state will be used in turn to calculinte 
physical quantities vhich will be compared with experimentally chserved 
quenti ties. 

Berne;les and Primers? (hereafter temea B.P.) have presented 
@ theory of solid te? ena have drmm extensive conclusions fron thelr ap- 
proximate ground-state wave funrtion. The present work is mmimiy an 
attempt to give a more ecourtbe Geocripticn of the ground state tnen the 
fommlation of B.P. 


Properties of Atowls and Solid Relive 
fold fie * poesesees no such startling; properties as the liquid 
phose. There has been no such intensive epplicetion of powerful theovet~ 
domi tools to the solid, as there hes been to the Iiqeld. However, much 
experimental work hem been performed, ond theoretical emmlanntion made of 
the observed results. Several revievs of the properties of aclid ne’ 
hee been publiahea.22495 





*c. Donb emi J. &. j folie eit we. Progress ar Low 





. Myeicc, Vol. IZ (C. J. Gorter, edito: Intersclence Publishers, 
# LOT Yee Sy e 


3H. Rernartes and il. Primekoff, Phys. itev. 115, 963 (1962), 
l | ; 
‘fF. London, Superfluids, Vol. 22 (Wiley and Gons, Hew York, 1994). 


ee a. Atkins, iiguid Hellun (Cambridge University Preas, 
Cembridge, Logo). 
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fanee 1943, He? has been made avediable in increasing enowtta, 
end there has been coneidereble recent work with this isotope. Rromd 


alecupaion of the properties of iiqtda and solid Ye” 46 contained ta 
& 43 56,7,5,9,10,22 


3 


Sve nuclear moment of Hey ie ~ 2.127% muclees mmgnetons. Ita 


i 
Me cole} ia 3.01603, while thet of fie ip L,.O03S0. 


atomic mans (on the 0 
Properties which ave directly due to the etomic electrons are the seme 
for the two isotopes agide fram an extremely mani effective mess correc- 
tion. These effects include atomic spectre and te interetamic potential. 
% eteospheric pressure, He, ike ie” » Tem@ins iiquid at aysolute 
nero tempelatisce. Crystallisation tehes phace on epplication of aout 


3 atm of pressure. A preasure of twenty-Tive atm is needed to solidity 


lie”. 


oe 
the properties of solid He” have not been thosoughly explored 
mperimentaily. Much nutleer reeonance erperlimentetion has heen carried 
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“a. 3 P. ions. _— —~ Derperabune wes of Meliun Three. 


Pr 
. Vol. I (ct. J. Golter, editor, Intem 


ty, 2 P. Peghkov and 4. 0. Cdinoveva, Geporte on Progress in Shysices 
Yolk. EXII, 904% (1959). 
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Liguid Heli: Three 






1997 


"Yoltun 3 
Colmbus, Chic, 
10 needs 3 0f the Sfeventi: ese gual Conference on Lov 
Temperatiune Shyeic: - 4, Greles . C. holtia Hallet, editors, 
uiveraity of 7 roto Press : 1962), p. “363, 


> | 
+". 8. Orilay and 8. F. Memmi, Ligaid end foltd He’, Progress 
in lone Demperotuve Paice, Vol. IIT (Cc. J. Gorter, editor, to te pub- 
ished). 





(J. G. Damt, editor, Chico Stete University Prees, 
Gi Leo). 
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out. WW. 4. Seirbens: bea reviewed this work. Dete on the FYT 


reintionsnips near the melting curve heve been published oy 0. a. Gridiy, 
ys i, 
. G. Sydoriak, ang 8. L. Mills,” and by Grilly and Miaie. Date on 


thermal conductivity heve been published Oy E. J. Walther ani i. A. fair~ 


Sehuch, Grilly end Mlis ia through X-ray diffrnection, heave Dcunt 
that selid te” eiets in ovo phases. The so-called a yhaae, 0 bexky~ 
senter cubic configuration, extets at preamues below about 100 aus, and 
the § yhese, oF hexagonal close-packal structive, exlats at higher 
prepeures. This phace teengition wag elso obsesved by Gritly end Mlis = 
es 6 systematic Ai econtinulty in the sywcific volume. 

Preliminary measvrements o: opectfic heat of the aclld have been 
nate vy Baveris, Boum, Brewer, Daunt end Melilliews.” he pressuze- 
tempexnture dependence of the molting curve han been obtedmed to 9.06" K 
by the smue group. they find the entropy to have o constant velue close 
to Rie over the range 0.07° K to 0.7° H. 

The principal enomaly in che behevior of the solid Mes in a wini- 
tm in the melting pressure curve, plotted ag ma function of temperature, 
et yp * 2.3 etm, T= .32° K. The existence of the minimm ic accompanied 
voy interesting effects, incinuding a neshtive themmnl eipansion cosfMctent, 
and the comsequent possibility of melting the solad by lowering ite 
terpernature et constant presmure. Ry the ceuat un Clapoyron oquati on 
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2 
toa. ¥, Schuh, B. a. Gilly, end oh. Lb. Mills, Pays. fev. 11), 
775 (1959). 
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14g Ba) 
the exietensce of higher thermumbmentc order in the liquid then in the 


polid wrid yield ouch behavior. 


she lieliun Intespartecie Potential mid 
piky Tem ltonion 





che stationnry states of an aggregate of beliua atoms are soluticua 


of the Schrodinger equation 
ee , B, Yn? (2) 


vhere The is the mea body Uemiitonian involving @11 the nuclei and ali 


the electrons of the eciid: 


K gt " H 4o Q 
me BEY ER he S 
20) "e as) ~ dyed “43 


eo A, ee 


Here 1 levels one of the © electrons, a labels cue of the 
H# nucleii, the X, Javel nl1 the coonlinates (Anciwiing opin) of 


electron i. WNuolemr coor@inetes are dewerived vy the y's. Y, and 


ek 
Ven ave the nucleous-electron and nuclevsemicleus potential energy. 
' . ' Ly 
Followiag the adiebetic approximate of Born and Oppenheimer, *” ve 


separate into a part dopenting ue electron coordinates that contains 





Ui Born end J. 2. Oppenheimer, Ann. Paysik, Sb, 857, (1-27). 





one 





nuclees coor? 


eoori nates 


nctes Ov paremeters, and a& pact depending oniy on nuclear 


BOW 6 ty My * -%,) Oy,» - 4 + Vg) 


We detise the part of tte depending on electron coomiinntes 


ona nuclear charges 46 ii Amen 
i ¥, 4 f “a 
e oy o + @ Voy (X, * 8 © Xo) “env, I Yay ~via Ty, Hye . “KL, 


(3) 
where wo show the Schrtdinger equation for the electrons, from vinich ve 


optain the electron eneriy and state, for a particular arrangement of the 
mucleii described by the values of the parameters (¥y> Vy} 
Through application of the adiabatic eyproximntion,”’ wheve ve 
include the fumetionai dependence of EB, omit y on the y'e, assuming 
only that the electron motion cen adjust at every ingtent to the changing 


positions of the nucleii, and with the motion still determined by (3), 


we obtain 
i >? 
. - G Ot Ely Vy + + + Hylh = {4) 
ee 


We see from thie development thet the internuclear potentiel EF . 
[for that i6 the role 4+ nlays in (4)] in a function of the poaitioas of 
all the misleii, end results from the presence of the atemic electrons 
and nuclear charges. 





& 
15, complete development of the edtebatic approximation applied 
to the present problem ia given by FP. Sekts in Modem; Theowy of Bolids 
(MeGrew M211 Book Co., Inc., Hew York, 1940) po. APO. 
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Do obtain an exapiieit faua for Eos detadued study of (4) da 
required. 

Wirechtelder, Curtiss and mn” nave deocribed the camonly ueed 
yhemomenoclogical fom of the interaction potential of He” atoms, and 
haye given @1 overview of the theory. netoes! hae Gomipared several of 
theae forms of the potential. more” bas meade a recent contribution te 
the theory of the repuisive core of the potential. 

At Gistences lasye compared’ to the Rokr redius, the interaction 
ef two heliws atoms as deceribed by (4) to charecterized by the von der 
Waals interaction, which arises from nttrnctiion between the electeic 
Upole mments mutually indueed in the electron charge distritutions. 
Tiis interaction energy depends inversely oi the sixth power of the sep- 
esation Gdietance xr of the etoms. Tiere ere also higher order tems; 
the next io the dipole-quadrupole energy, vhich is proportional to 25. 
At close distenees of approach, of the omter of two Bohr radii, the 
potential is steeply repulsive, mcoughly proportional to ~ COE (r in Jj). 
The electron clouds are stebl lived erovnd each nuclew: end strongly repel 
each other vhen they stert to overlen beceuse of both Comlowb regulsion 
end the Pauli principle. 

The phonopmenologioal forms of the potential have been chtained by 


fitting the theoretical Tosome for the potential to the tempecniuun depend 





ws. oe a gene | GO. F. Cuctios ma a, 3. Bird See 
theory of Gases and liquids (Wiley ond Sone, New York, 2958) yp b and 
A59f5 





“I3. DeBoer, Physica gh, 6 Ol, (2958). 
1s Moore, oe Cher. Phys. 33s y7h, (1960). 
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ence of the second virial coethiclient of the rave gases in the grseous 
state. A foo of the potential which hes been calculated with Low teenme- 
sture debe end virlch hee been found to predict the virial coefiicient with 
good accuracy to et leest 1° A is that of Thekingbem, Homiiton and Mooney 5” 
hereafter termenl 5.0.4. Taide potentiel has the form: 

“iGOr  h, 4 


re 
+ 


Vie} = | 770 e Ks _* erg G ©.61 A. 


(5) 
ah. 0 2.52 WI? | 
V(r) «| o77 © Or eS : | PE es eve oy & 2.61 A. 
P ug bike 


Anelytical application of this foum would be difficult. Nowever, the two- 
pareneter potential of Lenneni-Jones 


We) = be [(™ - BP] (6) 
With «a * 10.20° H, %2,.555 A, gives nh very close fit to the 3B.i.M. 


potential. its tehevior in the region of aevere repletion varies fran 
thet cof the more precise B.H.M. form, but 2% will be apparent inter thet 
this deviation takes plece in a region of such high vepuledve potential 
eo to leeve pryysical sesults umerectea. 
auxee beliun atoms in promimity undespo a three-wy Indurec dipole 


22 
interaction, which follows the relationship ~ 





We, A. Duchingheam, J. Memilton, and H. 8. ¥, Meese, Proc. Hoy. 
Yoo. (London) ALi, 303, (ioha), 


“Or. io inter tented thet the besic plymical effect of the 
remilsive core ia to reduce the probability density cbteined froma the many- 
poxdy wave function practically to sexo. Edther foum ef the potential is 
equally effective in this regex. 


aaa SWBON , De Chan. Piva. eas LOOT » (2953). 
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a ‘ ot RK Det + Bayes + ~ 
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where or 5 =e the Cistence between otoas a and 6 mengured in Bohs 


re 
redii, Boos 49 the induced triple Gipole exergy, and Baby 28 the tte 
body energy. In cases of interest to us ay > 2.28. Then 


Bene! (Ba, 7" ve + B,.) is of the onder a? >, and 69 can be ignored. 

Hence we ave correct in stating that Evy) 7. ¥y) ia @o- 
seribed adequately for cur purposes by (6), and the many-body Saniltonian, 
written to inclule all Coulomb and electronic effecta ta 





Y = ae a2 & 
ger -B ofed c vel a] . Cr |. 
ee) “ ~ @, 08d ¥ ny “ao 


Riise Memtitonieon is valid under the Toliowing liberi conilitions: 
1. Tempermiure lese than elxnit s6° K 
2. Density less than that corresporiing £6 a neares. 
nelghbas Giatanee of 2.5 A. 
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CHAPTER IZ 
THBORTES OF SOLID HELIUM 


fero-roint Motion in Tolid Heliuna 


We shali use the many-body Usuiltonien contealaing the suclear 
position coordinates 


7 ed L am f = 2 
with the potertial 4 


4 e | 
Megs) 





between peirs of heliva atoms. 

Cur purpose is to solve the Schrédinger equation for the ground 
state of the coliad, Yo and to obtain the ground state energy Ey > 
Given by 


+ 2 ot = 3 > . 
i Yo(X > K, ++ « k) z, VoiX,» Mp e+ s Koy) (3) 


vhere x, describes all the relevent coordinates of particle 1, inclu@ing 
spin, of my. Yo must heve the symmetry required by the statistica. 
3 | 


For He", which hae spin 1/2 end obeyo Fermi-Dires stetistics, we require 
that 


PY * {~ 1)” V5 
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waere P is an operator which permutes the indices of the coorfinetes, 
ant p ig the parity of thie permutation. The ctate in (4) t6 necessarily 
entisgymmetric under interchange of ail the coordinates of two particles. 
ite” obeys Bose statistics and has a wave Puwiction symmetric under inter- 


chenge cf the coordinates of two particles 
Dy «my 
Prove Yo * 6* AY Wy * Ue - 


We can use the ground state to ootein propesties of the salid at 
finite tenpereture by superimposing phonons, end opin waves, or opplying 
the moleculer field theory of permangnetion. 

The potential (2) t% cherecterized ty a week attraction at Large 
Be » @ minions near 1.1200, and o strong repulsion ineide r«=c. The 
two parameters coc and ¢ cen be assigned yoluea to give good vepresenta- 
tions of the potentials between rere pos atoms, and fair approximations o? 
other systems, ouch a6 nutieonenucloon, nitrogen-nitrogen, etc. Tenace ve 
ean consider « broad class of problems by reducing (1), (2) and (3) to 


dimensionless form. Following the Ciscuesion of Dernezsves ," we define 


Li, m By €y 
Ty &4 4/0 » (5) 
“~w aad om 
wk) @ ok o) = otF,) 
ed - es Pa F, 


EOP PRL LL MT ERAT, By A TOD OTD, RIA 


2 
Th. Rernerdes, Dyeory of Solid He’, Heliun Three (Sdited py 
J. Damt, Gyio Ctate Univeralty Press, Columbus, Ghic, 1060) p. 115. 
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The inat relationship is given by (i/o) ) O° /ak,” [Wy o)] 2 ae 2a a, “yla, Ds 


where the Laplacion is scivmeatically renmesented by the cecond pertial de- 


> asl ww me ow ed rw 
B= \ v(m, - * - y) * 2 aha, * @ - X,) Us, o 6 » Ry) 
= Pee 
= %, Hi, . . » Fp) (6) 
>. oe 
vhere AW Bo wt 
28 Gg 8 


Equation (6) aggiien to mny ayeten] syste detcribed by the two 
peremeter potential (¢). The physical properties cf a particular systen 
ere oll contained in the paremeter »% anti the definition of Ey: Velues 
of % Tor various subatences are: menon, 7 = 10>; neon, 7 x 10”; 
He’, 0.3473 ile, 0.302; neutzon-proton interm:tion ° I~ (20/(V.%,°)), 

(Lo ™ miclear radius in fermds, Yo = potemiial depth in Mev, »’ ~~ G.3). 

the ideal situation would pestit s genernl form of Y, with 2 
ag a continuous parameter, but no such ei:miified solution has been ob- 
tained for (G). 

In the limiting cage in which 7 << 1, where the qantim 


aspects of yy can be neglected, we obtain 


z, O° hE VG, re ee Cech (7) 


In this case, the enrzpy of the growl state is the potential 


energy lettice sum. The expression (7) is valtd for most solids, end 
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omiinary inttice dynmaics, inchwvling quentiged sound Slelds, stast wlth 
(7) and Newton's second Jex. 

Por Langer 4, the Kinetic energy tom of (5) (the sero-«point 
energy) ic an important part of the energy opevator. The theories of 
Uigquid helium end nuthear matter heave the common central uroblen of the 
handling of the minetie energy tema siuulteneously with the tvoeperticle 


potential. Solid Ne” ex te share this difsieulty. 


Sheories of Molid Helium Pour 

Several nuthore have comldernd the dynamice of solid te’. A 
vrie? deseription af their .worl. pertinent to our problem io given hore. 

(a) foownt® hes considered the shave of the potential around the 
lattice site. The tero-potnt motion emands the solid so that the nearest 
interparticle digtence lies outelde the miminem of the two-perticle poten- 
tial. The resultant lattice sum shows strong anharmonicity. Tooten ape 
promimnates the resulting well by a harmonic osclidator potential through 
& warietional calculation in the thermodynamic Sree energy. Eis results 
(ive the thesmodynemic properties of a veveroly enharmonic lattice with 
iazge sero point energy. Application of Iils remulte in workeble Porn is 
hindered by the necessity of maising Debye-type approximations. Cuantit 
syimetey effects cannot be ootmined fran hia model. 

(vo) eucher™ hes aioo considered the aetelled shape of the 





*} 
“D. J. Hooton, Phil. Mag. 46, 422 (1955); Phtl. Mag. 4G, 433 (1955); 
Phil. Mag. 46, 485 (2095). 


“" 
“I. J. Sucker, Proc. Phys. Gee. (London) 73, 965 (159). 
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potential well from the lattice sum. Ue expands the potential about the 
lattice point (0, 0, 0) as follow: 
2) 0) t 
We’, * Ws" ty te) 7 wi Cs: + ;" +2) + w(x F t o°), 
m ¢ 
thue ignoring terms Like sy" » 2n omter to heep the prcblen sopareble. 
In this expression, 


My wet, o(y) 
mm 
ves (H48(e) +E o*(2)] 
es te rid 


7a 0% er) + hfe F*4(2)) 


ig = 1anGo = (92) + G/e Me)I 


vnexre “(r) in the two-particle potential. She problem can be solved as 


G supexpooition of three one-dimensional anharmonic oscillators, using 


Wyo Va end W),« 
munker obtains an equation of atate 
Sean > 
OW Is , ? 


dHhiere Ey represents the seso-point energy, and is a complicated function 
of Kp, ext Wye fs Guoker points out, his veoult ie wmeatisfactory; it 
is besed on an Einsteia model mx therefore ignores the effects of netgh- 
vors other then through the idealised potential. ‘hic consideration is 
4n edition to the separability anprocimmtion made above. 

{c) Dugdale and Mactonsid* caleulate the onergy I(r) of a pair 





i 
‘J. S. Dugdale and D. K. C. MacDonald, Phil. Mag. 45, O21 (1054). 
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of adjacent atoms sepereted vy o dietance (4) ina rere gas solid. They 
Gomme eo Lennari-Jones type potential energy, with e first-order correction 
for anharmonicity effecte. They further aenuse thet the kinetic energy 
contribution to the total enermy ia proportional to the Debye temperature. 
Time Debye tecpereture £5 assumed, in turn, to be pooportional to 
(ae fax®)/ a. Using these assumptions, they obtain good predictéionn of 
pressure vs volume yvelationships and cohesive energy for the heavy rare 
goees, ond osder-of-negiitute results for ne”. 

‘hey find that 

Mom Vor! 3 


a - Le 14 x iD 3 
O 


Ms mole: mage, grams/rol 
V> * molar volume, om” frole 


B * molar exergy, cal {mole 


1/2 
es b/kGB| oy | —. 
eon wx nm ax K (Debye tampenature , Rs tattice 


pareneter ~ unity). 


The value 1.145 2 "3 matohes experimental values for Xe, hr, 
by 
A, Re end Be, with a mevdenm error of ot. 


” 2 i 
Tae wor’: of Hooton, Tucher,~” and Dugdale and MacDonald" shows 





CLleavly that thermal excitations cannot be treated by strenight agolication 
af the Debye theory of specific heats, om that the Inttice opecitic heat 
at low tempernturea muet come from a move fundementel considermiion of 
the dispersion igw for seund weves in the solid. ‘Their results inclule 


tie effect of the semp-point motion, out they do not include a description 
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a= lla 








eee se ee eee 
er rn ews ee ee ee 
et eS ee Om oe 
al a 


te us - 8 


ar es tae | 
im aii ~myn- 0 7 
Oe ee re 

sina <9 emma wong (62) a 

Nae 














oe et eter Ole ee Ws) eee 
- Mk a a Be 2a = ew 

ee a 
aa eke > Aue 28 © me we ee te Cee 
wim aaah Sa! Ee tiem 8 Ee OFS oe & 
*. (ere Lee 8 8 ee ee oe oe Ow oe 
an me mre ee 
eee (eee ee oe ewe ew eye &— © i 


~ Gu 


of more then one aton in the lattice dynemice (other than through collec- 
tive motion in sound weaves). In other words, they do not include the 
effact of ryytions of the nete@iibore of om particular atom on the potential 
of thet eaton that axices from the neighbors. 





meoies of Solid Heliua Three 

The theory of solid ile? hes been considered by Bernardes axl 
Prinekets? (hereefter termed B.P.) and by meesecten.” 

(a) Bevrnawles and Primakeff (Meltler-London approsc:) 

tn approaching the solution to (1), B.P. umed the potential 
(2), onl a variational weve timetion consisting ef a peodluct of Gaussian 
finctions centered at lnttice points: 
oF 


1, const e . o29,(ri). (8) 
id jo J 


ts §, 16 wsed to evnlunia “he energy expectation value, the 
sisynuesity in the potential fimetion yields on infinite result. To 
evotd this @fMaulty, b.P. mdify the potential by a factor wiich gives 


the limit sero es »Y, : gees to vero. Their effective potential is 


» \@ a \* oY 
Vogel?) # he En at rd ° (9) 





Ghe cut-off factor exple s a” te, re 20 chosen for computational 
convenience. The value ef . is entinated by calculating bnown properties 


+ 
of fe” using (9), with veetous values of c, witdl the best Mit in 
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ad 
“i. Bernaxtes ond il. Primekof?, Phys. Rev. 129, 968 (1960). 


On Demmerdies, Phyo. fev. 120, 1027 (1960). 
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obteined. In perticular, the cohesive energy and root-mean-square 
derivation of atoms frum lattice sites for tie wore calculated by a verin- 
tiongl calculation using the peremeter oa. The valve of cc thet gove the 
most consistent resulte in these calculetions wan then esoumed to mynly 
equally well to He’. Using thie esmueption, they cbtein the expectation 
value of (2) for the He? mae; 


K(¥, ae nf ar’ f(z") ek -| gr') 


Ti 
ee i ao ~, on 4 
+= A | fj @ a, 9 (z,) P"(e4) 
os am 
wet I ~ i, - 3) + rt) gr! ) d(r! 
@ lie HW) + Me U, (10) 


where z ic the coordinate of particle 1 in a froae with Ry » the 

position of the lattice point 2, os ovlgin, and the integration is per- 
formed over ali space. The specific voluze V enters through i, ~ B 5° 
The aomeption has beer made that the effact of Introducing symnactry into 
the wave function has a mami effect an Kk and ite variations. B.P. we 


the condition 


to find K. 
leving found a value for G as a function of volume, 5B.P. form 


& eymactrized weve fumrtion and calculete en exchange energy: 
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Jy. + Jy . (12) 


exe Rae is the vector between nemvest neighbor lattice points, Je: is 
the Kinetic energy past of the integral, and J, ia the potential energy 
part. 

B.P. find that J, = 3.2% 10 e end d= + 5.b 2 10%, with 
J= = 2.2 x Oe = = 0.08 K. 

Fron theiz values of B and J, 3B.P. then proceed to cbtain 
values for specific heat, entxopy, susceptibility, emi other themaxkaeaic 
properties of the solid. 

There are no satisfying supporting physical erguments for the 
effective potential used bY B.P., and for their assumption that the ob- 
served behavior of tte matches that of He? {through the conatent <¢ in 
the effective potential factor). The miltiphying factor in the potential 
energy exe” @ s™) 36 Gecigned to circumvent thie Gfliculty. Nowe 
ever, correlation effects owteide &.GA, elthough important, are not in- 
cluded. Other questions arise which are open to discussion: 

(1) ‘the product of single-particle orbitele used by B.P. does not 
include particle correlations peoperly. ‘This fenlt of the Gingle«~particle 


t 


function hes been diecussed by memy authors. ‘The difficulty can be wyier- 


stood by condidering 


AD vy 





2. % Heabet, Reva. Modern Phys. 33, 20 (1/61) Giscuszes the cor- 
reletion problien aml gives mew references. 








- 7 = 


‘ - - iS ht 
. e. tw. | . | i 
ee. e i | 
- ‘ cA —! 
Le 
il a A a we es ee 
- 

Ge FRA» ee POM | we ew 8 
wo wer- + Mew ear we 

ef oe 4 Oe ee oe 
ew Lee ce lee 
- Hae ots S colmaeng 
2 Sern owe Pee 
ne me cH i hn alleing STS 
ae —— ee ———— (ap 
Oe 
Phe oe cee et len ~~ 
oh Oe Ore ire Ge A me mY OO (oy 
eee) « a ae ee ee eee 

a RS of) Ree Cre eee coer er LL) 
milage Re a ae EE nn eed 


— ee ee KS oe ° FO |? — te 2 et 
ee ee ae 



































—_— Ce eee i" —z — a . « 
«= = pu oe be elo — 


CVip > * SG (2) Hl2) Vypleyp) G,(2) Gpledae, ax, - (22) 


if? the single-particle ocbitelis have eny region of overlep at all, ami 
Vio(7y>) has @ singularity atronger then | fx, then <V,. > is iatinite, 
contrary to Wiysicnl actunlity. Thlc consideration shede doubt pootieularly 
on the computed exchange integral, which Je determined ty the fPehavior of 
the vave function in regions of overlap. 

(2) B.P. performed theiz calculations essusing o face-centered 
eubic Lattice. They make the conjecture that theirs maerical valine ap- 
ply at lene’ omilquantitatively to the body-centerxl cubic contiguration. 
Rovever, since the exchange -integral is highly sensitive to particie sep 
eration, their coleculmtion of exchange is much Jess tthely to be appii- 
cable to both latiice structures at once then thety calculation cf 
cohesive energy. 

(3) The form of opin~inttice structure in an entdiferromngetic 
eryotal is densitive to the exchange integral with nenzest and mest near 
est neicqhbors 8 in & body~centered cubic Inttice, the next nearest 
nelghbora are only 15% facther removed than the nearest neighbors, and 
the exchange energy with both nearest and next nearest neighbors can be 
of the seme order of magritude. At temmeratuves high campared to the Neel 
texporatum, the susceptibility is given ty the Ciole-Weiss lew, 
xe of{? + TH). the Curie tenperatwe %, is given by the energy change 
ina opin flip, including interactions wlth all neighbors. (In Chapter ¥ 


there is further discumsion of megnetic effects.) 
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Oy, il, Ven Veeck. dd. phys. rei 22; 62 (2952). 
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(vy) Bernerdes (Nou-overlapping Wave Function oppronch) 
nervmvdes’?”” nes considered le” and lle’ from the point of 
vwlew of nonm-oveslapping sincle-particle orbitelea. The functions cly 
eve oquare-well ortitels, ani Dernardes calculates the integral 


vos eet SG) Pld) Vesey - Fel) G4) G,(adae, ae,» (43) 


USA ne 
1/2 win(x o 8, @} 
B.(t,) spay Teo Ee) e< a/c 
(Ae) 
G.(&) 80 §,.2 a/fo 


He obtains « serlen in (a/ac) for <ai> a <> + <> >» wad mini- 
mising <E ps » obtains a veive for nn. The series expansion coefTicients 
are aitferent far body-centered cubic and faxc-nentemd cubic Iattices. 
Hence Dernacties can deciaze one lnttice or the other most stabie. He Tinks 
thet both tie? end te? geese clonest paching, whereas experimentally tig? 
ie body-ceontered cubic, and ne* is benagonnas close-packed . 

fhe integral (13), elthough it ts finite, still overestimetes the 
potential energy. the true wave function gives a different probability 
Gensity describing particle “a” for each assumed position of particle “b”. 


She pimple product wave Tunction gives « distribution fimction for portichl 
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“NK. Bernasies, Phys. Rev. 112, 1994 (1959). 


Dy. Demmanies, Phys. ev. 190, let (1960). 
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"o", in a foma chowing no dependence on the posdtion of particle "b’. 
iis shortcoming pivevents an accwnte estinmte of average potential 


energy, except in the lintt of large seperation. 





CHAPIGn IiT 


A VAVG FUSCTION Wi CONELATICE: 
TUR PisNIRAG MOOD 


Meaning of Particle Correlations 
fn elementary discussicon will clarity the meening of the tem 
“yerticle correlations.” 
The wurnl description of a may-vody system in coordinate space 
is @ product of single-pesticle oxbitain (hereefter eniled 2.P.0.), 


a : 6%). (2) 


fal 
Here, the form of the $'s 1s determined by the boundary conditions imposed 
on the oystem and by the Schrodinger equetion. 

In (1), wo have ignored symmetry requiremerta. We mevely state 
that correlations resulting frou the Paul! exclusion principle are act 
ou: comcern in this aaa fhe product (1) in the veve function of the 
Testree approminntion end is the vave function of B. p.* (epert fron syne 
motrisationos). ‘This type of function can be used to deseribe exactly the 


state of Uo non-interacting perticlea, and cen be used as a variational 





uo be precise, this product must be constructed to give a definite 
Gyumetry under interchange of particle coordinates. The resuita of Shic 
section are unchanged if ve incivwde oo by writing $ 2 4 Pe, (2, ae 


where P pesrmutes the indices of i, » egsoigne the proper elke y an an) 
to the permmitation, agai normealises the product. The epin variable, if 
eny, is included in x, - the symbol oP, will be consistently used to 
repregent apace coormdiintes oniy. 


HW. Bernexdes and I. Primmkoff, Phys. fev. 119, 966 (1960). 
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eterting point for interacting particles. In general, (1) fedis for 
strongly interecting, dense syetons. 
lee, Nuang, and Yeog in their anelysle of stroaghy-interacting 
(hardephere) pesticles, we (2) in « particular forn: 
$ = ny itgpgenay MS vee By) UL (FL) Up(2,) --- UR). (2) 
Hexe the U.(F) are a complete orthonormal set of wave fimetions. The 
coetficients C carry the syometry of the wave function. Ghia Rmction 
om be used to calculate the hard-core problen, if many different conbina- 
tions of &, are used. Leg, Huang, an Yang transfer (2} to a represen- 
tation expresseé yy occupation magbers in i: space. They perfoon a 
verietional caiculetion, valid at iow densities. In a crystal order the 
@olution in k space must transform to = solution in r spece shoving 
eryotelline order. This operntion han not been successfully performed. 
HesiGes feiling to demonstrate cryctelline order, the theory cf lee, 
iveng, end Yang is valid for dilute syutems only and invelid for syste 
ee dence es the ome under consideration here. 
To help Qeecrlbe a ctete showlng crystel symmetry we keep the 
basic description in coordinate space. Ve use a less drastic separation 
¢ @ then (7) ond include e functional depeniense on the mitual sepera- 


tion of the perticies, 


« m $,,(8,,#,). (3) 
g ang oa’ 's 





3%, 5, leo, XK. Uueag, and C. W. Yang, Phyo. Mev. 106, 1135 (1957). 
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Ve study egumtion (3) by considering in sume detail the behevior of two 
interacting particles in e@ Dox. 

Tre usual interpretetion of the probability density corresponds 
to ea contour gap, with the frequency of obse:vation of a particle as the 
Glevation, ant the position of the observation as the location on the map. 
The “aap” should, of course, show conditioms in a 3-@imenaional solid, 
voativer then conditions on « piane. 

In Mig. l(a), we show, in a qumlitetive presentation, the probe- 
bility density for two perticles contained in so opherical well of infinite 
potential end radius v,. The curves show the “elevations along «@ Line 
azewm through the "map." ‘The state is represented by $,(7,) (8). 

The configuration fo intended to represent the steve of lovest energy 
(ground state) in the S.P.0. approidmation chbtained, for eocuple, by the 
Hartree-Fock method. If Fig. 1a) has validity, it must corsespond to 
weekly interncting petticles. We have clsmaky seen, in the case of 5.P., 
that, if we calculate ao potential energy velue useing 0 5.P.0. wave firme 
tion, Gfficulty can avise. in the integral 


evo es f OCR) 168, + F910 zs 088, 
“3 "2 


ae 


we ancune V = AIFS, . Por n>, the integral diverges if $, eae 
$ overlap anywhere. In the present case there ts comiete overlap, an 
n= 12, so that Mg lfa) cannot represent helium atoms. 

Ve vext improve the descciption by retaining 2 genesal Dimetional 
gependence on the positions of both particles, writing ¢ = w(t,; = 5) 8 


wie 


y(2, 2). We cannot sepernte (1, 2) inte a product of $,(2, 20,{F)- 
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Figure 1. Particle Correlations 
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In Pige. ib), 1{c), 1{a), and fe), weing thie function, we show hev 
the “contour map" concept must be modified. To drev the curves, ve must 
think of a vhole series of mmpa, whore each map corresponds to 6 nartic- 
wlar position of ome of the perticies, and displays a Gistribution 
function for the other perticle. (We think of the figures representing 
@ syste with the exponent n> 2, weakly attractive outelde 2, , © ge) 

In Mg. 1(d), we show the probability of finding particle ((2)) 
[hereafter lebeled simply ((1))J, at rj, [heseefter labeled (1')), re- 
gerdlesp of the position of ((2)). We employ the Direc 8 function to 
Locate ((21)) an ae, at i} , endurite fly(a, ay|" KF, “ #1 )a*s aAe 
= fly(2", 2)]°ars, = 2(2"). we plot ef"). 

In Pigs. Uc), Ka), and We), we illustrate the probebility dis- 
triwution of ((2)), wien ((1)) ie lecaliaed at (1'). This distribution 
ig gives by f{{y(2, 2)|° oF, - #1)a>(#) a ly(i', a)". The particle at 
(i') i¢ chown an a@ heavy vertical bar of hetat f{1'), illustrate 
the relative probability of each of the mans ifc), 1(4), end ie). The 
(iemeter chosen as eabsotasa is the Glemecter through (1°). Partiale ((2)) 
a8 distributed with cylintbical eyumetry svout thie dlometer. Hence, the 
novmalizations of Pige. l(c), 1(4), and i(e) are related in a complicated 

In conpteucting the figures we have aseuwd a strongly repulsive 
interaction of redtlus Zo" I? tae interactions were Wank everywhere, the 
representation 1(n) would introduce a omller correlation error. 

For thiwe and more perticies, we proceed in the seme femhion, but 


with a new set of mane for eeth additional particle. > ese le, with 
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three particles, Figs. i(c), Ua), and Me) can represent the conMigura- 
tion of ((2)) enti ({2)) for a given position (3') of ((3)}. We cae write 
the threesparticle fumetion as y{l, 2, 3), a theee-particle orbital; as 
(i, 2) 9(2, 3) ¥(4, 3), the twonperticle approximetions or an ${2) $f{2) 
$(3), the aingle-perticle appromimation. In the spirit of providing en 
improvement to the 4.P.0., and since we will have no three~particloe inter- 


eotions to deal with, ow work wes the tuo-particle epproxsimation 
O(2, 2, 3) = o(2, 2) of2, 3) ¥(2, 3), 


or, in general 


. . » 
@(i, 2, 39 « ® ® i) me tl WF,» uy) 


4 > 3 


Use of the Interparticle Coordinate 
to Provide Corveletion Bffexte 

fhe methods previcusiy applied to the ground state of solid ie” ; 
which include a description of the ground state (for exemple, 0.P.), 
express correlations through e product of single-particle orbitals. As 
descrived in Chepter II, such a description vYailsa in the case of a poton- 
tiel with & strong Singularity at the origin. Including the interparticle 
Cistence in the state function would remove this objection. The tuo- 
particle function of Chapter II offers a way of doing this. 

Previous work has considered this wroblen. in parallel calcuin- 


tions of the temperntiwwe dependence of the susceptibility liadd te? ; 
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& 
rue" ismtroduced a Mined interparticia distamee, aad Mee” exoured « 
hawt-core repulsion. 


Jestwou,” Aviles s¥ LWEEROCO Ch 


eel Deckman ” have ceveloped a 
veriationcl method (the Jestrow method) uning a trial Ancticn II 6, (x, ) 


Mf f (rz, ~ #,).~° Tia function resembles the I y(F, ¥) of the 
7s 3 4>43 


lget section. The Jastrow method includes a series expansion in clusters 


integrain similar to those used in ocotetning the partition function in 


debe 
clossionl otetistical mechanics of Genee svetens.” ‘The cluster amthod 


of obteining the veelationnl parameter ¢€ is wunulte?d va the salid stnte 





peceume 1% rasimes Iimiuid CABOTtery ond therefore the Jastrow method will 
not se:ve oux needs. lovever, the method has polnts of similarity with 
our method, axl we will retaum to it Inter. 

Berveckner’“ nas developed a method for incluiling two-perticle 


imteractiona correctly to arbitrary oiler. ‘She theory is a seattering 





My, 3. Peace, Phys. Rev. ST, 259 (1995). 
49. K. Mee, Phys. Rev. 97, 263 (2955). 
Ga. Jeatzow, Phys. Rev. 98, UT) (2955). 
13. B. Aviles, Jv., Aum. Phys. 5, 252 (1998). 
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Ome suthor is indebted to Professor Walter Kohn for pointing out 
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“5. F. Meyer and M. CG. Mayer, Statistical Mechanics (Wiley and 
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theory, end inclwies an expension similar to (2). A workeble wy of 
epplying the Beuecner method to s solid has not been developed. The 
paeucopatential method” Shaves this dloacvantege. 

Plawaés 7 in otudying the ntoaic mpecteun of helium, used the 
eet of coomlinntes introduced earlier by Hylieras. as Pluvinege added the 
feature of incliwiing the interperticie Gistance explicitiy in tho seroe~ 
order Hemiltonien. The sesc-order wave function Is thus a fuxtion of 
the space coorlinates of the porticles, conteine the intecparticle dis- 
tance, and cammot be sepurnted into « product of eingle-narticic Sunctions. 
it pooseanen none of the disadvantages aocribed to cther eethada in the 





Uaing Pluvinage's method, Waleh end Boru oe”? Chereatter toraed 
W.B.) hove considered the enerey levels of severel atoms and, in general, 
get better results in sero omier then the more complicated Pirst-omler 
ceiculations wruelly performed with hydrogenic wave fimetions. The method 
can algo te mplied to our problem, with remsonable expectetion of equally 


Good recuits. 


ihe Pluvinage Method 
Ane Pluviange method introthuces the interparticle distance into 


the meny-body Hanailtomtian and were function, ylelding an exmct foul 





135, Pluvinage, Am. Phys. 5, 245 (2950). 
Ly fA. er lLlerns » SMI GS Y nerpeee 65, oO") (21930), 


15 ' 
P. Waleh and C. Borowlts, Phys. fev. 115, 1206 (155¢); Phys. 
Rev. Lis ae) 1274 ( 3% 0), 
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exnprepoion which incluiles explicit Gependence om all the interparticle 
distances, as well en the spetiel positions of 11 the partich 
Ve write the meyy-body lami ltonien: 


Y ws h | 4 a 4 


Mere we Geseribe 4 identical perticles of mass m, interecting with 

the potentiat V(r,,). We have already econ that we cen take Yr, ,) to 
“VJ 

the Lennend-Jones potential 


Vir, ) = dal Cafe) ~ (o/s)? 
¢ # 10.00° £ 
o* 0.556 A (5) 


Viz) ie en externally applied potential, the seme Mmetion of space 
fox all pexticies. 

We will now follow Phuvinnge in introducing the interparticle dis- 
tance into the Hamiitonion. 

In any calsnuistéen involving n independent variebles (x, 2 Vy» 
Wy Use Yas Woo s++), We oan introduce auother set of m variables (a> 
Yy Bis Aye Vn? ey cee}, rOleted to the rs oy "7 a (u'), aad poouibliy 
olimplifty the calculation. We o@n even heve m >n, ena ey 
alno provide m-n 4Znienentent constretnts in the functions < u af 
In fgpendix A we shov that the following statement 16 time: 

x 9(0) sz y,(a,, Vas Wye hye oo We)» mn fumetion of the Coxtestan 


coorhinates (uy, + Vio Y } of NH perticles in a cyotem, we may eet up 
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2 
another function ¥p(2") * ¥ol%y Yy> By +++ Fy Syr Myo Fag ote Mya? 
which conteina (8)(N - 1)/2 extea coomtingtes, and we mey require that 
4,2) ss ¥=") everywhere in space. fn order to do this, we must have 
(3g + (8)(NW ~ 13/2) independent velations * s (ur) » end the rule of 


aifferentiation 


Os  - : det Mp 


pup Ch owernpen 


Su uy dw su! se 


In our cage, the damit itenten comteins the interparticle distances 
M44? O06 We ate led to introdkme them as entra coomliinetes. 
We then have ‘ 
e353 4” Ss SU 
(6) 


i Rh. Etec oa a 





«odie ee a 


| 5 a on a 


monte 


an the 3M + WC - 13/2 independent relations. 

in what follows, we ohell .we ec weve finction of the tyne > " 
which givea the game remilts as Vie Simply because we require ¥5 a3 ¥en 
We stall find thet Yo yielde G& wievul neturel separation. 

Ve drop the dictinction betveen iy end Vy » but remeber the 
apecteal role the x's play in the finctionnl form a? F4o° Fron Appendix 
A we heve the rule of di ferentiation: 


> 3 al 
e“& * 2 ee 3,5 ?) 


Vere, to the rignt cf the arvov, we ppeak of a yy, typo of dmeorlpticn, 
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Pron Appendix A, we heve 


”~ 
“ 7% x Sy 3) 
a . ee ‘ f re y 
EB A anh 4 7) * Bt y +} we (- ees 5 "i | 
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«> o> 
Qe. ., V a > 
+p at Oc BR (8) 
if} "4 44 aft!’ ~ ie "4G 
Vere, we have yar 
“ © s  % a 
V4 See & * 4* eo &£ (9) 
scien * ae "e243 
& x ~ x e 19 
8.43 (x, 7, (19) 


_ oA => 
Here 84 is the vector % © Oy It is used to express relmtive dis- 
placements in tems of the variebles x, © U,, t, “V,s 2, “ v,. The 
symbol 2, 4 le one of the IM -~ 1)/2 aefAitionel coordinates, 
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ce ae nN 
Mas “ Mr, 2, ) + (v, “ v,) + (v, rm w, ) . We need doth oe and 3, 
ms 
Since derivetiven in x spece, Like 3fsx or Jy 444 2 Operate on Suse 
wo’ 
CZ0NB Of My Fy By Ondo Dut aot on functions of a y" idinewilse, Gerive 
atives such as ofr, r Go not operate on functions of =, but do operate 
on LunctZons of v, 5° in mldition to the internerticle vector B, 4? We 


will also lave need of the position of the canter of muss of pais if, 


expressed in XM, ¥; % Gpace. 
on wh os 
2%, 98, +5, (41) 


We oan write the Hermilitoniean UH, 


Heh, + a (32) 
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An attempt to we the total Nemiltonian tc detemine « set of 
eigenstates of the fchredinge: equation eppears wiinviting. Intend, ve 
we an epproximation introduced uy Pluvinege. We lock upon HU’ as a 
perturbation, and look for the states of Ep: 

i, $ se gy (15) 
Bs Bas wit Fy) 
Pips Pyge ts Fags Mone see Pays oer a We atteck (215) by the 


Here § is a fimcticn of the 3% + MM - 1)/2 variables ( 


stendard method of seperntion of variables: the form of ly makes the 
solution separable as follows: 
G= cel OF) a, | . xa) (36) 
P 4 . 2 © Os L>d a3 
Here we show the epin function ye © fumetion of the spin variable 
t.: the state fumetion @ i9 en approximation to the true-stete fime- 
tion. If the perturbetion potential ' is truly small, the approsime- 


tion wiil give tittle exror in the energy level of the ground stete. 
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B= Eh, + <B i219 > (27) 
<O€1G > 


In working with the fumetion §@ , we must remenver the role of 
paremeters of integretion. We must, therefore, follov oa progran which 
has the following steps: 

(a) Solve (15) for § 

(o) In uweing BR, ond it', let 37x operate on «x, gy t, 
only, and let Pr, F operate oa =, j only. 

(c) After differentiation [step (b)], introduce the constvaint 
yy Uy J + (yy = 2s? + (a, ~ 

(a) After introducing the constraint, evaluate matrix elements 
by integration. 

Yote that the omer of steps (b), {c), and (ad) cennot ve altered. Aa in 
eny ennlyois izrvolving a paraneter, we can, hovever, Cirst introduce the 
constrnint, and then perform any remaining diffoerentietian by ordinary 

rules of differentietion not involving 4/9r,,- 5; «(We shall find that this 
latter prograxa will permit the calicwintion of the average cohesive onerpgy. 

Pluvinage end W.B. heve shown thet (17) ylelds excellent valves 
of © for several systems they hove analysed. Purther, they have guc- 
ceeded in obtmining good resulte by considering only by 3 etter mering 


edditional acsumptions comverming the function *. 
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Renovel of the Singulnr Interparticle Potentiat 

The Gemiltovien (11) etill conteine the singuler interparticle 
potential. The various applications of the Piuvinage method differ in 
the handling of the potential v(r, j ). Here, we shoy how to uae the 
Pluvinage methad with etrongly repulsive interactions. 


ib in (15) contains the following terms depending on r, 17 


2 
ae 1 > 
sp sj-Bv® +dyice.)| om xe,,) POR). 


We can make this tem vaiich, end vith it, the potential Vi g'% a5) by 


muuietng”” 


wp 3. | 

+ Bi “r,ay * 2 ve,)| “A,’ * > — 
Expression (15) can be locked upon es a mathematical artifice de- 

signed to sumsoumt the difficulty of the aingular potential. it om. alno 

ve considered a Schrodinger equation deNning a particulor member sof the 

comtinuum of otetesa, the Lowest lying of such states, thet of sero energy. 
In Appendix B we solve (13) wy chenging from variable r to 


variable q«@ c/s, teking ’ as ae pover series in eecending positive powers 





Wie could write the right-hand sige of (18) as ea4 Ye g'Pag? 
We would find no bound states, hence Sy zO. For G54 > O, there 
would be oscillations in 4% at Large distences, which would cawe real 
Srovubie in desexribing on extensive solid. Since we hope to have the fumc- 
tional fom of “(r, ,) be independent of Ty49 we munt Olininate oacli- 

oa 

letions et the outset by having ¢, g * O. We could otild solve (10) with 
CG, ie j on the mght, and then set ¢ ~ 0, but this process would oniy 
serve to introduce additional tems in the series expaneion for “ wich 
would contein ¢€ as a miltiplier, a needless camplication. 
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of q. The reoult, for the plysice. peremeters of fe’, is a rapidly 
converging sevies. The two linearly independent solutions are 
“9 * Fyf é/r) 
t, * C,(9/r) B4(o/z) 
B(x) 2e gente KG - 7+ b> 3) x + BO 5) x” 
Oo ~" has 


oO + 22-8-7-%- 30a Fe... 


Px) = 1+ fb. 5) x +a 8+ 5 +b) xo + BM - 20) es... 


o> 
Here Ban de K ifn, whith hes the velue 15.0 for ie”, 


At lege rr, x Cos aid when rw«9o from above, % tends 
rapidiy to sero. Thats behavior is cheracteriatic of a haw! cove interm- 
tion. At large r, hy oO; it then reaches a maximum st an intermediate 
value of yr, and heo @ fintte value at rec. In Appendix B, yy is 

shove to yield a solution with an outgoing Tix: at large distences. ‘This 


is en unsatiofectory characteristic vhich cannot be permitted in $,- ve 


thus choose Cy = O, 

We now denonstrete that Cy mx i. C2nce the Hantitowiean does not 
contain the time, we may choose X to be vonl. fince F is ven, Co 
must be real. 

i 
Cy 23 Co : 


At lasge Gsteances of separation the interaction haa no effect on 


the twoeperticie wave function. Therefore, the probability density for 
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particles i and 43 is in@ependent of r 44 ei multiplying the veve 
el | 


function by the square of the correlation amplitude leaves the wove func- 
tion wichanged. Tiat is, 


a 2 
hor ~» 0) Ae Bs, ~co)} m i, 
"ys " ry 
Therefore, Co Fy (0) #1, or co * lL. We teze the tive root of 
2 
Cy ® 


If we were to wee the WAD epproxnimmtion to obtain the veve func- 
tion for r<6o, we would find thet %* would continue to fell off 
rapidly as 

fe 12 


__ 2 . 
exp| - f [Via - v)2a] bh” a@r| x [¥(o ~ van] . 
reo 


Hence we are justified in esswaing that <X(r)} is venishiugly omell inside 
r © o (provided only thet we use coution in calculating derivatives for 

Y & o). Ovr representation for 4, the series F,( c/r), has poor con- 
vergence properties for r<o. Because of the anliness of <(r) for 


r<o, we con cut off Fy(o/r) by writing 
3 ee i; - + 
Ax, FY (2, ¥44) alr, , &%)) (35) 


Here ty 5 c where x, cennot be much less then o, end @(x) is the 
atep function @{x)= 1, x>0; O(x) 0, x <0. 
om 1 
2 a “Cr 3 r 
In Fig. 2, Xia ys x (, 5), wnt =X ag) * (ry ts 5 are plotted 
Lor’ te? - In Fig. 3 we iliustrete the cormerigon between the semies Pinc- 


2 4 
tion Fo and two approximete engiytical representations. The function 
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Sis is the trial fumetion vhich hes been found to give the lowest energy, 


ont of meny tried in the Jestrov neti. 


@naova in both Piga. © and 35 te an extremely close Tit. 
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CHAPIN: IV 
ME MANY-BODY PROBLEM IN CCLID IBLIUM THADE 


Density Metrix FPommileation of the Mony-Docy Wave Punctica 

In this chapter we shell develop 2 theory of the many-boay problem 
which can be applied to dense solids. The develomment uses the density 
matrix formmeliom of Lowiin hs which ve veview briefly here. 

We coneider the meny-body wave fonction yy which ?ulfilis the 


symmetry condition. 
P yf a4) a oo 8 8 i.) Sd ( “3 _ W(X, » 465 > & «© oy) (1) 


vhere P is a permitetion operator vorking on the subscripts of the co- 
Ordinates, ani p i6 ite parity. For fooe perticles, » *=1. For Fermi 
particles, ) @ +l. ‘The spin variable, if any, is included in E. A 
physical quantity & enssociated with the system ion represented by a 
Homiitien operetor ©. . It may be expressed as: 


op 
o.°@ #24 42, Ba. B Oe... (2) 


wnere we sepactete the operator into ite ger, one~, two-, . . . particle 
operators. The prime omits terms having two or move equal indices. 

ZO evaluate the aevernge vniue of 2 op , wo introduce e@ series of 
censity matrices of various orders: 








ib. 0. lowain, Phys. Rev. 97, Ur7h (1955); Phos. Rev. iT, UO 
(1955)3 Paya. Rev. 97, 1909 (1955). 
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ei =, 
(Xj {%,) = mjy"(2'23 . . . WH) y(ie3.. . . M) at, . (3) 
The symbol ax, oigaisies integration over ali varlebles excent x, - 


rk! xs ix, x %) S ME > 2) f o¥(u’2a'3 . . 2) e(ies. . . we) 





ck, a, 
@é e e ps e i Lf @ S x = 
(x! % m5 3 Key » 


ie * tr; ¢ $ : 5 ” or v7 
sae pr fee eee fF cee it) (iz eer JF seve B)(ax, ves * 
ft 
(23) =s aw >,;2 - ve ae ’ +) 
P(X] Kye e Kile, K+ + + Ky) ey (2'2'3" 22.) ylaeg ... HM). 
The netrices (4) Sulft11 the following conditions: 


MK, Kylky Hy) = CR, RlRK) 


and 

Nx! xk, %) = 9 NK x ad bo , &). 
Dat {Ky KNX, x )* Mk! X! RX, % Ky)» vagerdless of 1. (5) 
Ve also have 


(pe), 2s = tt & - 
r (% i e @¢ @ nly Jo . @¢ @ Bed! 


bed — p) ye et 3 a iz - — —_ 
Hsirepl (xX ees rd, x13 ty Shy eee 5 1 3%, (6) 


Sne dlagone? elements are of special importance 
Cyernhle 
¥(X,) @ y{x,]2,) 


%,). (7) 
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Snese matrices are symetsic for 5 221. Ke con interpret r(&, av, 

as N times the probebility of finding a particle in dv, at Fy, with 
opin ov regardless of the positions end apine of all the other par 
ticles. Sintlarly, MX,, %)av, dv, ie the muber of pairs, MH - 1/2, 
times the probability of finding one particle in av, at ¥, with spin 
g, ané mother in dv, at Tp with opin §, regerdlese of the spins 


and positions of all other particles. Accomting to (3), wr hove 


frye, = as fra, eek, &, = ww - 1/2 
rx,» A, “6. « Xai, ai. ees <x - pili - pe (8) 


> ; 
The matriu P) is entiaymmetric in each set of indices when 4 2 ~ 1. 


Therefore, if ta cr more indices ave eqint, If P) #QO. In particuler 
—_ —_— 2a —_ — 
MX, %,) a, 3K, , Xs R,) a 0, (9) 


These relntionsivipa are expressions of the Teuli principle. 
We now consider the use of operators to obtain average velucc of 


physical. quentities. Using a two-perticie operator as an exemple, 


: "in o Ww at, ... &. 
Evaluating the intepeal for the generni peir 1, ve heve 
<apeean-~1syfo,ym,...& 
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Agplying (3), wo get 
ca>=f fy, nz} a, He stent Sead am, m&X, (40) 


1”? 


We use the convention that 25 operates only on the woprimed coordinctes 
R, end x 4? end not on x: . a ence the non-(iLagome] clomente are 
needed to expredc the expectation values of differential operntern. Using 
the convention of (10), we get, for the general operntor (2), 


ss 


afe«* a ; 
< Mp7 S¥ o,,¥ GH, . + . Hi, 


t 


= *- — > 
“a Ja, r(Xj 15, az, + f 2,,NX 
Sh ice (21) 


in our problem, we heve the two-pecticle operator ¥, l=, ° ®,| dy 
end do not Introduce theee-perticle operptors. Hence, a description of 
the #0130 in terme of density matrices of omles 2 should give accwsge 


values of <2 "Ve as follows: 


(i * x: 
Ln >= J (% + Le] a, * +(3 3) 0) MR, lk, Ret aK, (12) 
oP i NR, ek, c 
In (22), we teke care of nowmmligetion by the dencinntor, which is con- 


structed to give weity for 28 on al. 


nangle-Perticle Density 
(a) In tite section, we we the fommlion of the density mtrices 
divtxexhied in the preceding section tc build up an NH particle descrip. 


tion of the solid. Giese deseription will be anpprostiaste, in thet it will 
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deseribe the properties of the solid in the single-particle gense, just 
as the Giater detezminant deseribes weakly interacting sywtens. Our de- 
ecripbion takes the correlations explicitly into accowl, an additional 
refinement not incliwiled in the Slater description. 

(>) Gingle-Perticle Probebility Density. ve shall write dam 
the f-body density metrix in tems of the fw Phovinage weve function, 
end reduce it to a single-particle density metel«x. We then show, for 
Gense wyotems, that a simple ropresentation can be used ag the solution 
to an integral equation, yielding an explicit form for the single-particle 
density matrix, which includes correlation effects. 


ane U-body density matrix is 


(Ht) . - P 2 
r lz P 4 $4.) . By Mey) (13) 


in (13), we show the Glegonnl fom and do not introduce peine coortlinates. 
fhe result is scurficientiy general for our needs, since we do not sea: 
expectation velues of differmiial operators, but merely the probability 
Gennity. Hon-diegonal complications will concer us in inter sectiou. 
When ve use (13), we must romenbver the set of constreints Ms a 8, sy 
it io not necessary to aopecify the form of “ until wo attempt to make 
inmerical calculations. We shall ascwae, hovever, thet “~ is a universal 
fumetion of the disteace between amy tvo particles, independent of the 
opin of the particles. 

Proa Eq. (3) the single-particle denodty matrix v(3, ) eon be 
ontelned by imtegrating the f-paerticie metrix over all. coomilinates oxcept 
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lil yen 


r(x,) = (const) TE... Ks. Ra, (2b) 


The constant is a nommaliaing constant Gepending onijy on the number of 
particles. In the first section of thie chapter, wo evaluated it es N 
iteeir. 

We cm rewrite (14) in terms of the represestation (13): 


r(2) = (const) © 209(2) $5(2) S452) 6,(2) -. Om $00) 


a) =" “ ‘ ie 


a "3 = 
%53\"a9 “onan? © ° + Mae By Gy Bs, oy Byes - (15) 


In (15), we have emplicitly ghown the diegonml fom of the density matrix. 
The pemutetions refer to interchanging the coordinates of the G. Sox 
both ¢; and $, cimultansouniy. Tie resulting function is syanmmtric in 
the interchange of ony two particles. The symbol P.. (5 for symmetric) 
ohowe this syometry. ‘the lest relationship in (5) shows this property. 
ia asmmption is implied in writing (15) in thic wey in that we 


are neglecting the opins of the particles. We heave assumed that if 
' gC") G82] 10) 4,0) 

2) (9194439) ara | 
G32") Mle = 16.2) Olet, 


r®)(aa}ae) = o.(at? 16,(erl® + fecorl® 1 Car’. (26) 


By eo doing, we heave dropped fron M2) (a0} a9) m1 overlap tort: 


- [(oe(2) 9,2) g,(2) %,,( 1)) + (complex conjugete)]. This term: cormt=lbutes 


" 
positive oz: negative enounte to ~ 2) » deepening on tise seletive spin 
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functions of $ «ona ¢,.- We ahead find that (a) correlations keep our 
particles well nepavated, and (b>) the b'n are orthogonel for particles 
weil saperated, so that tributions of overlap are very anall compared 
to terms Mike 19, (a1. Nence we ase justified in cropping the cvorlap 
term when it appears in (15). We wlll find thet we camot deo this tem 
when we consider exchange effects. 

We ahell aeowee throughout that the syametioic permutation of 
@eronal elements shown in (45) is ademurte to doncribe probability dens- 
ity effects. We ghnil find later that exchange otlecta are negligibic, 
thus justifying the aovengeice. 

Wratang (15) tn moditied fom, we have 


v(2) = (const) 2 [1B Py Oo°l2) 0. O:°C9) oo GAD 
5 -™ 


The tem in curly brackets is piBeh); 4 , sow OF dg 5 * Ay vos Bi: 
Because the diagonal fom of © is oymaetric on intescthange of tt6 ca- 
ordinates, the omlering of the labels 1c imantecial in the $'s cs well 
en the X's. lence, in each member of the maa over §, the multiple 
integrais are identical in form, independent of the position of the index 


j. We can, therefore, factor owt the inteprel in each om over J: 


we fe & 4-1) om "? 
r(2) (const) ( 2 4,3 y) sx ($5 2 ee OO My Be ON 


We must now consider the central assumption of thin development. 


He aggume that we can project out a single-nerticle description from tae 
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angy~particle description (13), ox thet there exists ea yx{2) vhich de 
scribes the average behavior of the solid try ansigning certain properties 
to ao single average particle Imbeled (1}. ‘The vetinement thet we build 
into y(1) shoulda include the effect of correlations between the particle 
With eoowlinave ry exi the average position of all the other particles, 
Dut will not include correlations within other pairs. (these other patra 
would enter into a two-body density matrin description M1, @), which 
describes (a), correlations of 1 ané 2 with the averege positions of o11 
the athe: yperticles, avi (b), the actual correlation between 1 and 2.) 
The Gul hap) in (27) casry the correlation offect between etan 
1 and the other atoms k, wat we need to mow the everage positions of 


{ y.. 
the other atu. For this purtose, expresaing pi ) 


ap a product of 
gingle-pesticle density mxtrices 1s adequate. 


Wy 
u x 
x Ma, @, 3+... B) a mt v4), 
| jan 
we have an expreasion vitch satiefies the relations (3), (4) and (5) aw 
plying to density mateices, and specifies the positions of the varticles 
wlequately for the purpose of evaluating the integral in (17). 


We con now write for (27) 


a oy = 
2} = pomes 3 OO7(21Y nm vit) tf Mile) & 
r(2) : ®, (2) a" ) 3, apis.) &, 


(36) 


» ij 


iy 


¢. 


2 : 2 
. aaa e, (a) A J Us) Me cK, 
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Bguation (30) ie an integral equstion in y({j) wich gives the 
wrojection of @ ningle-perticie atete from the moxry-body etete. Fo ab- 
tain it, we have eapguned: 

(3) She two-particle separation of the weve Sumction, (13), is 
an adequate description of the No body niysten; 

(2) Mmong the single-particle sintes, overlep has a negligible 
effect an cinghe-particle properties, in that integrels of the tyve 
f£$(2) (hav, & # J elt ventohs 

(3) X 3h? 4) ig @ universal Ametion of the distance x, 4 
segemilesn of the identity of i and j. 

We nove not sasud ro comliitions for the 6's, nor a 
epecific Punctional foma fou %, Aon 3° 

In the nent section, we snl) use (15) to obtain en emlicit 


single-particle representation for solid Le. 


A Golution of the Density Equation for Dense Solids 

Tt seeas obvious that ¥(,) in (24) is e function which hes a 
meoximim in ite spetiel dependence at lattice pointe and shows the transais- 
tional. symmetry of the lattice. Ih fact, we can define a solid ly requir- 
ing thet it be described by oa r(?,) with these properties. We mate the 
eiditional essumptionn that to a good ovder of approvimetion v(F,) hos 
omerical syymetry in ite behavior in the Immediate neighborhood of Inte 
tice points (this considerntion is independent of the fact thet the 
moxie of the y(P 3) exe Cisposed according to lattice symmmtzy), thet 
viz) decrenses monotonteally mwyy from inttice sites, ax thet it is 
continuous with comtinuouws derivatives everywhere. A simple faction with 


these properties is 
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a= .- & 
“ (Rt) 


ze ° 
d 
— — ; 
Here By i@ the position of the lntijce waite j, end =», is the poi- 


tion of pewticie 4. 


We can dettine a quantity y (&, ), ‘hich is @ component of y(r,) 








ORY 3 
v,(%,) a @ -* | 2 
r(z,) ‘ r (2) (19) 


Ve rust remeber that (2, )° hes mewring only neer lnttice sites, which 
neang that Large values of oh, ~ 3) ive vs z,) "OQ.  Benee, (29) 
Presents no complication from overlap, since ve have se >> a » ond 
the only tem of y thet contributes near R, 2.8 Y;.° ’ 

We redefine tems in (28): 


r= FICO £4045). (20) 


The eubseript 1 appears on y,(1) since (20) only describes the density 
= —_ ao 

neex &B,. In (20 ; ip the distance 2, * * frog the Meld point 
4° In Js ty, the e a, 4] fra the Meld poln 


—_ 
- 


7 to the lattice site 3, and 


£,(u,,) 8S rF,) (hla de, - (22) 


J 
Putting the acaumed solution (19) inte (21), we have 


e(R.=,)° 
o yo o 
£j(4,,) = (const) ze  * J 6x 


ie “« ) ie 7 ow 
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Noting thet with O° > > , all the clesents of the sun over 4 are 


very eamll exept for t¢ = j, wer get 


= 4 
f(y.) « (const) f e “ (me, 


eo’ os’ oN 
where ds an y Y.s 8 coorMlinete with ovigin at the lettice site j. 
v 
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We partially imposed the constraint on r,, in evaluating (43). 
We now complete the process by evaluating G. for ,. ma ®, » the posi- 
tion of particle 1, remembering thes ? 4h 4 3) 49 a fumetion of the 
Giatence from lattice site j. The evelvation of G,. proceeds by per- 
foxming the sum ; 4, end r e, Zor constant Ay » corresponding to 
@ chell of atoms at constant radius, and adding over shells. 
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Figure 4. 


lattice Sum for lwo-Dimensional Square Lattice 
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We Gia not have to assume e Draveto lattice or spherical symmetry 


elintnate the tesn in cos @ leading to the result (29) 
Ve now have from (27) 
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prenaion for the probability density of on stom, determined by correintions 
with ail the other atoms of the golid. This important result is the pro- 
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Go wovk with (34), ve must comsider the 6's in more detail. We 
eooume that they fit the syuewetry of the lattice. For example, ine 
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where n@H4+3,8, 7-1, ...6,;, 5, 4) 3, 2, 1, 2/2, 2/9, 1/4, .- . 
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Evidently the veves with n<i can be used to conttruct a packet show 
ing the probebility density near ind‘tviduai imttice sites. The long waven, 
with n>i, Gs not enter Into such muvhets, but appear ven we consider: 
ceveral, particles together. To: exemple, ve can see from the dagran 
that if we consider 4 etemp in a row, we audt deel with n= 1, 9 = 2, 
n® 3, end nh. Longer wavelengths will not enter. 

In perfortting the integrel (22), we could just as well integrate 
over Ek end evaluate the correintion vith 3, as perfor the integral 
over a> 3° thie conpideretion, togethers with our forantion of the product 


we (ty 4) leads w to the following interpretation of af rn), 
d hs | 2 wel » 


a) 
It forms a packet of waves, with wavelengths shorter than about 1/3" , 
FS 
none Tt, Therefore, the remaining ¢$,, which appear explicitly, are 


Geteruined by other considerations then correlations. For example, they 





nh fe ere ee ee |, oe ee ee 
i pee a UE ee ee Le Oe 
hes oe ee oe eee ee ee 
ee ee ee ee 

ft OL Om om tee ee & ot © 8 eee | eee & 1 ee 
— we Ae ae Low 

eee Ge eat ee ee ae a? 
ke 6 ae 
——— <5 ee ce eee Oe ee 
+ PPP bo oe es a SO +l ahls 
gg We ee ee ee geen — > ers eo 
ee 
wot yee tee ee es ee 





cowld be called upon to permit constructing states conforming to the 
Dealt peinciple, or to mike up eigenstotes of parity. 

Lowoin ; in the third peper of his series on many-peeticle aye 
tems, describes the wee of two such etetes, those for ne 1 mad 
ne#i/f, as © oymwtrised trial Auction for considering pais interns. 
tions in a ayntenm such @6 gUuxD. 
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vhere of ia @ tem of onler ai/ 3, nese negligible fos: lage oyotems. 
In Appendix D, we dexive enprescions Lor G We Pind thar Gy, 
hae two parts, the fiswt pert, Bo» % independent of 8°, ana 
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Ghe part of on which Gepenmis on 0% is 








ett (Nil tee oe Ades eee 
ener it se 


re rmen ge ire tae 
hes ey PPT cars se — = = = a 
mt © et oy Oe Me eee 3 Be 


cere @ eo Dee = cum 
“te oe ee yg © eee 








60 - 


Dae 


1.0 


\ — £ from series representa- 
OFS: dO 


On, Of Na 


\ me 25 from Yukawa representa- 


‘ tion. ot x 
0.2 : 


0.02 
0,07 


0.005 


0.002 


0.001 ft —__+——— 
3.0 h.O 5.0 6.0 7.0 8.0 


Figure 5. Elements of the logarithmic Sum, Eo? Plotted Against ladial 
Distance. 





why Leo 


gE 2 
cae (ib, = r ny ~_* 
a9 le | 2) 
+ J#R + $OIn, - 23) adie 


where $ io the ervor integral. ‘there are af@itional negligible tems 
fy 
af order (ty 3)" > vxich are Giscussed in Appendix 5D. 


Pron (35), (37) and (39), we cen calewate 8° © By 


+ oa *s dy) 
By ME Moye * By) (mi) 


5 ® ~ ho 


We fairet obtain 0 by evaluating © os @® quantity indepeadent 
a 


. 2 Mm S14 2 
o? by » We then cvelucte . on te witi2 a soneistent vebue cf 5,, 


ined. Figuwe (6) shown the result of auch calevintions. ‘The 





48 abt 
ceiculeted values fell on a amocth curve which can ve fitted, with no 
discernible deviation, to the formule 
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Pormule (42) ta cherecteristic of Ho? in the body-centered cubie inttice, 
end im valid in the vange 3.80 <a < 3.00. Appendix TD to a sample cal- 
qulmtiion for a point on Mg. G. 
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in thie section ve uwe the single pamiicle prooability density 
we heve obtained in the last section to determine the average energy. 
The expressicm contains the Inttice constant a through the probability 
density permeter 6. 

Using operetoxe introdawed in Ug. (12), together with the deasity 
mxtedx of (13), we can write the W particle energy density matrix in 
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useful do differentiation, as explained in Appendix A. Ghe pevoctetion 
cperntor P' permites the prime coordinetes accowling to Pemil stetistice, 
just as © pemmites the uperined cooréinntes. 

We simplity (43) uy noting that the opermtor 1, ,* 1/2 V4 4 neken 
no combxibuiion, since (f, , * ife Vv. 4) % sO for alli, 3. For ow 
symmetry we ecounm that the tem in if' is negligible compared to the 
tem in @ As expleined more fully in Appendix F, Mrsetly H’ ics short 
range because of the Gerivetive aPr, which veniches et large r, 4 
and therefore ' 4s proportional at most to 7 end not MN - 1)/2. 
Hecoxtly, in a Dravais lattice, a lnrge part of T' emcelo out because 
i’ 498 linemr in relative Qi splacesent » end a Given atom has a povitive 
displacement with respect to one aten, and an equel and opposite dicplace- 
went with reepect to mother aton. A contribwtion to BH’ remains vWaich 
is proportional to the everege ogquare of the deviation fran the Inttice 
site. In Appendix we prezent a pleusibliity argument that this con- 
txipution 16 considerably leos than < f 7. 
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r(2) = 216, (2)|? 2 2,0u,,) 
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We must coneider most carefully the way in which we combine the 


conteibutions of the sai teatiai” ee ng on cach of the 7 y? When we 
employ the representetioan oe @ . ‘this special attention was niso needed 


* 
in eveluating £2 A, oad 4S in the preceding section. 
J > 
A, for 7" 
ib, oC or 
“1d 5 
to ect woon. Tach of the terms in the integm. in (50) introduces o 


Bach of the / :, hag a single variable 


“kinetic enexgy” of relative motion obtained vy taiing 3 ® hn,” ‘ (motion 
nomanl to the Line between 1 and J gives no contrition), and hence 
represents one degree of freedom. But F not only tekes derivatives in 
the direction Ry, » out also two more derivatives of equal yniue in the 
ta directions nome2 to Ry 4) or would trent each member 4} as Ghee 
pres of freeiom. Hence, to use the ophorically mymmotric matrix 
- 8° an” tO veprecent the influence of correlationn on the enerpy o? 

wus 1, in the single particle approximetion, we must take 1/3 the result 


obtained by use of the lanplacian: 
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P ~~ + 0 . 2 a= 
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+ J 5 OFA) 4,02) exp(~ ay” 0 /2) = expl- a,” 85°/2) af, 
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Ae fee x f6,(2)1° ae. (St) 
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In discussing Bq. (49), we heve already seen that the first tem 
of (53) 48 of order oF evi hence of ovder a 6 » ond negligibie if 
we neglect overlap aml exchange. 

Therefore, the single-particle energy, from the second ter of 


(53), 19 





Tice ez ae 
<B> = ge 5 5 a i oy " (55) 
Exchenge Pnergy 


(a) In this -c clon, ueing the Pluvinage form of the wave func- 
tion, we concider the sywmetry effects on the magnetic propertions of the 
solid. 

Intirely analogous with the treatment of clectron euchenge, we 
oagume that the Paull principle opplies, so thet we are limited to a 
state deacription of any poir of particles which is antlaymmetric in the 
interchange of ali the cocriinates of the pair. We attempt to find the 
wave Punction of two-particle stetes of lowest energy, one symmetric in 


the space coosdinates, correspomding to the triplet spin stata and the 
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other entigysmetric, corresponding to the singlet spin state. the energies 


of these states are conpered, and we obtein the exchange enercy. 


jy ead Bay > toa 


ret 2s the singlet spin state energy, any is the triplet energy. Posi- 
tive J cozvesponds to peraliel opine. 


Our caloulations will obtain < E> (yn = +2, -1) end a nommolize- 
2 
tion integral for each case, from which we can construct 
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Zhen wo heve: 
< BOD 
Cle = A ) 
z* Sng (55) 
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<i> -sF5 


Imchange in the mony-body case is somewhat different than in the 
two-particle case, for which we have just defined J. We define enchange 
ao euchange with respect to e class of neighbors, for example, exchenge 
energy with next nearest neighbors. 

We will compere the energy of the solid in tao cesed-«one, vihen 
the spins of 1 enc the olews of neigybore are parallel, and, two, vhen 
these epins are antiparnilel. If the exchenge energy with reapect to a 
Given set of neighbors ics positive, that exhenge integral 1s tenmed 
ferromagmetic. 

Ye proceed by weiting down the Nemtitonten. We delete the tern 


[1/2 v, ,* 24] because it yields sero in the product with %,,. Accord- 
ted 
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ing to Appendix P, i’, 4 will eppem: egain an a symretelcolly omtented 
gum, with the linear velocity denendence giving sero to lovest order. We 
take vV, (x, ) to be a constant. ‘The Nemiltonten (for Lowest omler exchenge 


enexyy) is then 


sar 
° \* t% te 


Mere, * includes ali of a class of neighbors, like, samy, nem 


sereat neighbors of the amme spin, which emmmot be distinguished within 
a 

our oyiimetry of space and spin. We have = 
x 
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Ve now get up the energy-denelty matrix, as in the preceding 


section 
Van(2 ) Pani? ) Vay i8 J+ 
i' @' oy, 8 = 
. 2 og” d,, ( ) %,, ¢ 9, 3 
ro ob wall a, + 2% a » | 
i wa . $. (2°) & (2°) @ (3') = =} Ieome, unprime 
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G@LlL the other particles. We have no interest in the symaotry of, say, 
2 end 3 as @ pair, since ve do not propose to distu'd their sayrmetiy. 


We use a form for $, here which bes o naturel symmetry for ox- 














change : 
WX rh A th ae 
ro) a _o- ee . a= 
$, (t) ~~ Ce CUS ae o0| ‘ cos| a] (59) 


each veve hes o phaae to give the grestent symaet:7, vith no nodes nt 
lattice polnts. Equetion (57) is written for e cimple cubic lattice. 
For other lattices, similar functices can be constructed. For exaaple, 
in a Natl type lattice, where the I veprese its one spin system and Ci 
the other, ve would have to gets of equation (%)), each with W/e nen 
pers, end each with a + (2//3)ea, the length of the cube edge. 

We cen rewrite (50), using « Sleter detexainant which conteins 
otly information about the pairs 2 - n in antioyometsic states, but 
which contains no temas giving antiaymmetry to the patra j-n, (j #4, 
ni). In rewrlting (50) we mst inclufe the other atoms in c form ate} 
quate to construct the contribution to the provebllity density of the 
oolid fraa the petra jg -n. Ve cen do thio by including then as 


3 19, (2)1° (4 42). We oon use the Slater determinant formmt, vutting 
at 
sero for nll the tenas not involving overlap with atom 1. For example, 
for taree particles, | 
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vVheve A» (2) ¢,,(2) (3) 
(2) O02) 6.03) 
c= > (3) (2) (2). 
if we edopt the convention that we drop tems showing overiap of three 
atoms, the product simplifies to 


AA+DB+OC-AB- AC + BDA- CA 
= [2 ny Iag(eal® WCort 14c30l? 


~ (9%(3) $.(3) [O02) (2) O(2) (2) 


+ )(2) O (2) $02) 6,(2)] + (aoe, 2 +3, 3 2)). 


Hence, we can use the determin. t P to consider exchange of 2 


axl 2, or 1 and 3, but not 2 and 3, nor exchange involving three 


QLCIE « 
We, therefore, write 
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if we understand that we include ovevlep between 1 and only one other 
etom at a time, for each member of the class A of ind@istinguishabic 
menibers, (GO) gives the same result. Using one of the matrix elemmrts 


out of A to vepresent ali in A, we get 


Orci) = af2") 8") OL") -- 
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We calculate the evermge valua of the operntors using « different function 
fox eoch, hence ve have a Gifverent normalicing constant for cach. 

We see that @, ds going to coutribute <7, 7 to the ouwa, the 
value obtained for the direct enezmgy. We can evelvate (61) for av velue 
of A formally, and obtein a formin. We oan apply the formula, ome 
Sor each A, @dding the reoultes to get the total energy chance in change 
ing the apin orientation of 1. 
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in (62), we heave used only that part of the determinant wich descxlbes 
the ayumetrized wave function of the pair (2A). 

We now adept the convention thet the index identifying the moin 
of @, }, 326 plus one for the spin of a opposite to A, and mime oe 
for the spin of ea pernllel to A. Zhen, the expression 

EO (2) OCA) % 9 OCA) 02) (63) 
is entisyeretric in interchange of 1 ond A, ani tee Poult principle 
4c obeyed 2f we we (63) in conjumetion with any symmetric fumetion of 
1 am A. 


We con vewrlte (G2),°reprosenting all atoms, except 1 and 4, 


by © single-particle denatty, as we Gid in (19) awd (45): 


4+ &, oo y 4 PEGECA) GAD +g BOD gFa)) i»? 
heh 2 Hy V Sal 4 


= OD) O00) #9 O09 6,(2)1 5s Fags) Gy a) 


x mes’ :  % oy. (Bay, ) cK, ), (Gh) 
with 
5, * Ota, a * Bisa) (65) 


She factor 2 in (G5) arlees because of the convention thet the 
exchange enexgy ic the Gitference in onerry for two particles upon exchenge. 
woerens (64) given the difference for one particle. Formule (Gh) ie the 
remult we would have obtained if we hed agewned thet if ato 1 ant atoa 2 


give om oxchenge energy J 1 ten the total amchange ener for 22 
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neighbors of type 2 is A Sans AD ny in the mudber of type © atone. 
Hovever, under this simple esmermtion, the question arises: hov oan wo 
fnelude the icinetic energy of etan 1 simultancowly in all the pairs? 
The derivation (64) encvers tris question by showing ws that the idnetic 
energy of eton 1 does not affect exchange at alt. 

ihe evnlvetion of (64) in ae major task. We con woe harmonic oscil- 
lator vave functions for the @'s. ‘They wlll inclufe the cut-off effect 
of the product over k «et lerge Gletences from the peir (2A), and will 
pemlt the syadetry to be affected by 4, if we teke one of the states, 
say A, %to be the first excited state, and the other, soy 3 tc be the 
ground state, in a variable desoslbing the degree of freedom clong tae 


an 


line joining the lettice sites gf the peir. he fimetion * can be 


“ue © 4 
A), A= 19h, b = O.Oob/A". Ze we 


represented as * = (l-Ae 
hanvitle ¢@ and ~ in the way decerited, they yiell elementary Caussien 
integrals (with complicationa introduced wy the cut off in x). ntro- 
ducing the behavier of the predict over ik in the region betuwen the 
tice sites introduces o crucial problem, in iteel? wortivy of a separate 
otudy. 
nstead, we calculate an wocimatve J, valid ot low densities, 
ucing the single-particle density matrix. 
Ye Jearned in the last section that we could got <7 > vy putting 
y{i'/1) = Jet") efi) v{l) (provided only dhot At + oe” Cann /5..) ole 
ommJ2). We obell use this iden to pedwe the density matrix in (64) to 


the Gingle particle form. We cen interprot the veve Ametion 
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4(2) 4,001 4 (2) aT vt 7] vk) &Uik) Ok) ak to. - be the emplitute 


for particle 1 to be near lattice gite a with particle A near lattice 


site A. 


Thais ie because (a) the product over kk keepa the perticles 


neer these two lattice sites, and (b) (ry, ) prevents then from being 


om the seae lattice site. We recall thet e repregentetion for (1) vas 
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=f a R R,-#,)° 


We want to bulld up & noneGingonal two-particle density 


matrix from this representation, end still have a resuit vhich shows the 


necessary conditions: 


mita'|i2)-e r(2'a'je1) = ~ ri'a'jex) 

mia}i2) = y(1) y(2). (6) 
We receil that ve had 

r(i) * 7 ¥4Q) 


8° (R Lak A 
¥ (4) = (const) e : 


The enticymattric form of the density matrix which we seek is: 
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[oeme, prime -» uprime}. (67) 


Formuix, (67) yielda the con@itions (66), and giver the correct 


Gingle-particle effects. 
“4 
we wece to use - n° /(on) vy for 


ao in Bq. (52). 


Pormula (67) would over-estimate the energy 1: 
the enexpy operator, in the same way 


We must introduce the same factor of 1/3 here. 


We shall use the form in (67) es a firet approximation valid only 


Go long as the syhevically aymmetric approxion 
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The spherically symmetric anprostimation is 


results cen be expected to have 
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Here b=af/2. The s tude joins the two lattice sites 1 amd 2, end 


&, end &, 


tice sites (that is, there are Inttice sites at 2° ; 


identifying @ with A” in (6h), 


the form (G4), we got 
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We get, performing the operations: 
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Yience, in Bq. (55), we hove 
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Ve we Ba. (56), but Include the factor 2 fran (65): 
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Beglesting fe > eompared to <1, we beve 





gs ~ 3 lat ew -e +eje ° ™ 
2 2 
a 2 a” a - . 
7-7 ° nx i (69) 


Equation (6¢) cemnot be expected to be valid in the high density 


case, ince the overlem of the functions yf) enters strongly in (69). 
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& better appvomimation would comsies of waing lese then spwrical 
eymaetry in ¥,00), foo exemple, 
om " = £ ¢ 22 2 2 » 
+6, [(x" + 7 + 2”) * Bo (ey eee” + xe” )] 
¥,(1) = - ; 
lowever, 1f we enticipnte that these corrections to the exchange 
integral ave of the seme sign as (65), then (G) gives e lower Limit on 
d. 





CHASE V 
COMPAGTOC TH EXPO ST 


Blastic Sivecta 


in Chaptec IV, we have odteined an eopveestion for the cohesive 


energy of solid te? in the oody-ceatered cubic pheeec: 


w) 


a . 
Bere & > (a) 


with 


‘- « 1.27 © &.Ga(e « b.092)° a (2) 


She expression for b." is obteined by constructing the grounl-stete wave 
Section for the given lattice structure from the correlation function 
Ar} = of o/r}, calevlating a," for various e, snd then Citting o 
curve to a quadratic equntion in a. 

Formula (1) gives the ccheoive euergy directly. 

Goodkind and Fairbank, , ex etch” hove observed the tanpernture 
dependente of the nuclear regonence relenetion times. They use the 
Disesbergen. Cimsell en rc ua theory of the diffusion relaxetion mech- 
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“d. M. Goodind and W. M. Fairbank, Heltwa Three (Daited wy J. 6. 
taunt, Chic State University Prese, Columbus, Giio, 1569) p. %. 
- A. Mich, Helium 7 (Baited try J. G. Dawt, chio “tate 
Universi 7 Press, Coluaous, ry 2960) p. 63. 
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anion. The relaxation times Collow the Arrhenius lav, T « 40 g " P 
where z, ie the activation emergy of diffusson. In Fig. 7, we have 


plotted the experimental date for the activation energy of diffusion, tu- 
gether with the cohesive energy given by our theory. 

So a SMret approximetion, the cohesive energy is the energy to 
create a vacmrry. Once a vacancy exists at iow density, little energy 
is required for the vacancy to move between “earest-nelghbor sites, be- 
ceugze the new’) nearest neighbors do not intrude into the space in which 
diffusion takes plece. “Te er, the aiffusion coefficient (en! -2e relaxe- 
tion times), chowl’ we determined, in the main, by the number of vecancies 
through the Doltanenn velativnship 8 = Hl, exp - Boao / (kT). Neve EL. 
is the erergy to create a vacaaty in o dilute system of vacancies. Then, 
if Des, D=), eup(E/(«t). The c.sve indicates general agreement 
between Ue diffusion dyte and our theoretical values of cohesive energy, 
3f we esoume tot c, 23 Bec’ At high densities the diffusion active- 
tion exerg-ea Lie weis sbove our curve, incicating thot our mechanisen of 
Giffusion should include © density-cdependent activation energy. Cur 
theory cannot treet this problem. 

Ve have attempted to prefiic. the trensition to ea hemegouel close- 
pecked strucvure at high deapities, using our spherically syemetric single- 
particle Geisity. This in done Yy calculating Whe cohesive energy in the 
two structures, the structure with .oweet enengy Gensity being the more 
Steble. Ga regulte are equally valid for hemagonel close-pm).e. and 
face-centered cubic in the syherically oyemmetric approiiwatio. Tits is 


@ better approximation for the 43-fold syumetry of the face-centereé cubic 
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lattice (with respect to menvest mi next nearent neighbor) than cor the 
Gefold symmetry of the hexagonal lattice. We Cini the clome-pack 

tices to be leea gtabdle than the body-centered cubic for all densities 
comtarncy to experimental findings of a transition to hexagonal close- 
Pached structure at high density. Uowever, the spherically symsciric ap- 
prcuimation fatis to dlaccininate betsmen hexnegoual close-pecked mul 
fece-conterad cubilia lattices, and it is therefore necessary to put the 
ayommetry of the lattice intc tae single-particle density in omer to 
meke & theoretical study of the transition. Because of the expected 
Similarity of close-packed cubic end sphesical syametry, we find the Tace~ 
centered cubic Imttice to be exe luded et a1) densities, compared to bodiy- 
ventveread cubic. 

fhe pressure is given by 





Pease. 


Hance , pa ob = 


centered cubic Anttice. Therefare, 
a | 


But Ve 0° /1.30 3 Jeter for a body- 





in Figure 5 we have plotted pressure, vamed on (1), (2) and 4}. 
fiso plotted are iov~temperature experimental pointa reposted by Eihards, 
i , ry 
Bawa, Brewer, Dewrh ond McWilliams and Cyferiek, Mills end Grilly.? 


~~ 





Ry. O. Raweis, J. i. Bam, 0. F. Bewver, 7. C. Deirt, and $. 6. 


MeWilliens, Ueliun Three (Edited vy J. G. Dewst, Chio State Universality 
Prese, Columbus, Ghio, TG) >. ms. 


; *3. G. Mydoriek, R. L. Milde, ena B. h. Grilly, hye. fev. levter 
en (1960). 
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Bernarites and Primazore” (B.P.) do not give a theoretical value 
of the prescure. However, the tabulation of cohesive enexyy shown 1n 
Table Ii of their paper gives an abrupt change in pressure Sram positive 
to negative at om lmttice constant of about 2.7 A, vith ea negetive presa- 
gure et high densities. Cuch behavior is not consistent with the somali 
emowit of experimenter data. 

We con calculate the Debye temperate, &, throug the bulk 
compressibility using en approximetion employed by Beranrdes. ' 

abe bulk modulus is related to the elastic constmits 


ev » : (c44 + 2 C15) /3 (5) 


We cannot determine C1 end Sao separately; however, we cer acaine 
that we have isotropy (cay -%, 22 Cra)» and that the Cauwly velation 
holds (e505 = Cyn) We then have 2, * agl3 ° (cy fo), Were - is 
the density anc oc t and ¢ q ‘See the long tudiaal and treusverse weve 
aout velonities. We eleo have f= 5 c,,/J. The Debye temperature 6 


ia Gefined uy 
eehe (64° af)/3, (6) 
. 
where 3/2" we fe* + 2fe,. Cae therefore gets c = 0.64 ce 


we now combine the foregoing assumptions of Deras des with our 
value of E to calculate 6. We have 





ited 


6, Bernardes and H. Primakoff, Phyc. fiev. 119, 960 (1961). 
tn. Bernardes , Lye. Mey. 120, OO7 {1560}. 
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5 a 3 > 223 Q i. 3 
But =p * o(1.3/e") BA = 2.3/e", BV e = “Sa 
—lUCmUmUCmMSE 
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pa +58 “3 We. 2 el. 
03% 30° a &. 383 
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> > - 
a = 4.86 j.2 = « oA *~ (e@ in A). (7) 


at? 


en 
Since E= 0.0) + 4.02@(a ~ 8,064) i A 











@ = 4.06 Vice fa - 2 (a = bo6n)}a - A 


@2 22 % (9 im a). (8) 
va 


There are no mxperimental date giving the specific heat. Uowever, 
near the melting curve, we predict a Debye temperature of 2).2° K, under 
the aseusptions concerning the treneverse modes which we described evove. 

We can aloo check our results with the findings of Dugdale and 
Mac Donald 8 that the Lindemann melting formula gives the some vaiue 4 
for all inert gases 

2.8 3 
uo V.- 


ne 





21.14 x 10° ig) 


Fa he Ck 


“as S. Dugdele and D. K. ©. MecDoneld, Phil. Meg. &5, 611 (22:94) 
(footnote 4, Chapter IT). 
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Mo @ moler cass, grems/oole 
Vo @ molar volume, on fale 
h, ™ molar energy, cels/mole 


@ © Debye temperature, ° K. 


Ve can eniculate a value of A Yor very low temernturesc. We 
tome a @ 3.0 A. ing our value of _[ and @ for 9 = 3.0 A, ve 
obtain 

A= 2.6 » 19°. 


this value is about twice that for other rare-geo solids. Nouw- 
ever, a factor of 1.4 in 6 coula cause the discrepancy. further, He? 


has extreme quantum effecta, and (4) may not hold. 


Dachenge Effects 
In Chapter IV we have obteined an expreseion for the exchenge 
energy of aolid te? in the body-centered cubic configuration, expected 


to be valid at iower density. 





@ .& 2 
= "* S' a 2 
jt. & 5 C . (30) 





che experimental data of Meyer, Adams and Feivbank” can be used 


to estimate the Curie tempereture Te as a function of lettice constant. 


Their date are shown in Mg. 9. We estimmte TT, by using the Gefinition 


a 
i 
or te contained in the Cusile-Weiss lanes” 








Saleen ae 
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“E. Dwight Adena, H. Meyer, and W. M. Poirvenk, Heliua Three 
(Edited by J. G. Daunt, Chio fitate University Press, Coluwtbus, Chio, 1960) 
p-. 57. The author io indebted to then for making their data available. 


19 + : r - ae ad 
J. K. Ven Vieck, J. plys. radium 12, 262 (2951). 
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Figure 9 


SUSCEPTIBILITY OF SOLID HE? FOR VARIOUS PRESSURES VS I/T 








(12) 





where % is the magnetic bulk guaceptibility, C is the Curie constant 
for te, 7 io the temereture, and 6 49 the Curie temperature. ‘the 


Curie temperature io given by the exchange eaercy 


Reve J, ic the exchange exergy {a negetive number) with respect to the 
cless of neighbors A, and nm, in the number of A-type noighbors. 

We heve estimated be 
mental ante.” We obtain the values by replotting the susceptibility as 


for various values of a from the enxper.- 


1/X egeinct T. The low-pressure curves thus plotted have en umeistak- 
able straight line portion, with a negotive intercept on the @T axis. 
This intercept givee the valve of ve » provided that the tenpereture 
of the portion of the curve from which the extrapolation vas made lies 
well above Bae In Pig. 10 we show euch & plot for c presaure of 35.7 
atm. In Fig. i1 we show the results of eeveral euch plots from the dete 
of Fig. 9. The extrapoletion process introduce. . uxerteioly which ve 
heave indiceted on the Mgure. Qn the seme figure, wo have plotted oe 
obtained fron ow theowy. Zhe eereement with experiment ia gocd at inrge 
values of the lattice constant, but obviously breaks dova for denser 
solids. ince cur theoretical value of J ignores the non-syherical 
symmetry of the locations of nearest and next nearest neighbors, and dees 
not include the effect of next nearest neighbors on exchange, one would 


ne led to introdkwe cubic aymmetry into (1) insteed of spherical 
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a 
ae Extrapolated from ate of 
Adams, Meyer and Tairbenk. 
a 606.0 atm 
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lattice Constant, a(A) 


Figure ll. Curie ‘Temperature vs lattice Constant 
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oyermetry. There vould then be a hope of predicting the scherp brem: in 
Z delow 2 * 3.6 A eg @ result of introducing interactions with next 
nearest ueighbors. Calculation of J, useing Bq. (IV-G4), would also 
predict the contribution of next nearest neighbors. 

We cen construct a pictuve to expinin the experimental drta, as 
fSollown: 

(a) AG lozge separations, correlations are relaxed aad there is 
ccpasntively lerge overlep of the atorstic weve functicna. ‘Tho exchange 
imtegral is characteristically negative becawe the lattice structure 
end correletions tend to hold the atoms apext in a region where 
ave, AX, , is positive. The symmetric space configuretion peraite 
the atoms to be closer together, with a lower potential energy. The 
kinetic enexzgy is changed but litle between the ctates of differing 
epmnetry. (Hote thet here we do not rely on the weve function X* to 
carry the intermetion, but refer to the interparticle potential.) 

(hb) As the atoms are brought closer together, correlations be- 
come “tighter,” the region of overlap between nearest neighbors becanes 
ameii, and the exchenge integral iikevise becames omalier. 

(a) At otill Lower velues of Inttice constant, the correlations 

tween nearest neighbors become so “stiff” as to greatly reduce overlap. 
lowever, corréietions betwoen next nearest neighbors are not as reatric- 
tive, end appreciable overlap (and alco exchange energy) will exist 
between then. 

All the preceding appear to be substantiated by the behevior of 


our predicted J with chenging lattice constant. We can extend the 
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picture without violating what we have learned from the theory. It is 


we yl. for antiterromgmetic subetences, with decreaging temperature, 


CCN 
to shov a temperature dependent increase of the quantity “T over the 
constant value of the Curie lew, Y= Cc. The quantity %“T veeches a 
maxime at the Heel tenperature +" At iower temperatures, “ST Yalis 


19 sal. 
+4” won 


to servo linearly. Thia bebevior is predicted by the theory 
interactions with next nearest nelghbors are taken into acsoumt, in a 
structurms where alternating adjacent crystal planes carry identical epins 
tn @ach pinane, but opposing spine from plene to pirne. 

Further, Lt is known wet, when the sublattices retein their 
mutwel cormelations (e.g., nearest neighbors evtiperallel, next nenreat. 
neighbors parallel) that “(f= 0) = 2/3 K(f = T.)- In this exprescion, 


T. dia the Neel tempernture. If there ig no conmection between sublat- 


i] 
tives (o.g., nearest neighbors without correletion, next nearest neighbors 
antiparallel, ox ordering of the second kind), this veintionehip dees not 


hoid. 





De. 3. Govter, Proceedings of che Qyrposiin on Soild and Lguta 
| te Un 


ree (Chic State University end Air Force Office of Sclentim@e 
» 1357) pe 1. 
12 

‘The Weel tecmerature D, ie the temperstura below witich an 
enti ferromagnetic gubetance dhova magnetic omtlering in vero external mey- 
netic Pleld. The Curle temperature Te i6 defined as the paremeter Te 
in the Curie-Weiss lew, Bg. (11). In general, Sy < Spe with their ratio 
(never greater than 3) determined by the relative magnituder of the ec- 
change interaction between nearect neighbors comered to that bet 1 next 
nenrent neighbore. Ven Vleck! ena 2. Tegteisn, X. Youwlén, ana 2. Kubo, 
ee in Phypies 4, 1 (153) give extensive reviews of eyed Servceang 
nection. 
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Reference to Fig. 5 shows that the date for $3.0 Atm have a 
G@etinite maxinua at »% = (const) x 5.6, T=0.1° X. The expected vaiue 
of X(T =O) 4f T° O.1° K to then X* = 9.75 (const). ‘The observed 
walue et 0.0625° Kis % = (const) 5.5, and falling slowly. Hence, the 
date at 69.0 Atm might correspond to an entiferromagnet with o Néel 
temperature of 0.1° K, with interactions vetween nearest neighbors and 
next nearest neighbors of equal importance, conforming to the theoretical 
prediction” that X(P = 0) = 2/3 X(T 2 f.)- 

At higher pressures the overlap between nearest neighbors vill be 
reduced still further, and interactions with next nearest nelghbora 
should predominate. ‘hen the relationship (f= 0) = 2/3 X(T =) will 
no longer hold. ‘The data show no tendency to follow this relntionship 
at high pressures, supporting the conjecture that the exchange interaction 
with one of the classes of neighbors is very werk compared to that with 
the othe2 clase of nelghbors. Our exchange cnergy becomes small at high 
density because of ommil overlop, resulting from the strong emeclusgion of 
the corr@éiation function Ne et emall Gistences. It can be expected 
therefore to represent interactions between nearest neipiibors. Tince it 
continues to fail as the density increases, we may state thet the theory 
indicates thet J (nearest neighbors) is small at iticgh density. Zien the 
observed marked increase in Tf 


C 
chenge interaction with next nearest neighbors. A magnetic transition 


at high density should stem from an ex- 


might occur when the nest nearest exchange integral grows large. A epecific 
heat peek at the trensition temernsture, which might be a comparetively 


hign temporature at elevated pressure, could explain the difficulty in 
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attaining therunl equilibrium reported by Adana, Moyor, and Pai rom: o 


Conclusions 


Using 0 two-particle seperation of the veve Smetion, we have 
meade an approximate solution for the ground ctate of solid helium thaee. 
We use a deacription of the solid in a oingie-particle renresentation 
which 1s determined ty the many-body properties of the solid. Ve ade- 
quately predict the density behavior of the cohesive energy of the solid 
in the body-centered cubic phase. Near the melting curve, ve predict an 
extiferromagnetic Curle tempernture of the omer of 0.01° K, and at a 
slightly higher denwity, a Curle temperature of ebout 0.005° K, in agree- 


ment with experiment. 
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APPENDIX A 


TNS DERIVATIVES oF TUR win FUNeTION 
WICH UTEPanvIcn Panera” 


We use the following, well-known theoren. 
i 2s") al f(x, 2 vy? ad Ry Yes? ¥39 " @ * ert won? a "32? 
8 peat Sf 2 : 
Migr 6 + + Foge Boys ts “neue ig go (mctéias” of g verlobles, 


4 


end if the x” are each Gitierentiadle avwilom of mM vErieoee u, 


i | ‘ * 
x = “> Va Wy Uy Vey) Tess @tce, Tres 


. re s 


au! 402 Sa ra’ 


in our ease, we have a spetial wave fimction gle’) Giu,» Vie 
Wye Uys se + Vis Vy) s where wu, V, ond wv are the Cartesian coordinates. 
The number m is therefore 3N. “he riuvinege wave function is a uev 
function £(x) = Hay) Fyr Bye Bye © + + Tyr Bye Pap Myge © + + Mog Bay 
pes He2.,8” We moke ench x o fimctdon of the o 





Ohne worl of Pluvineage and Waish and Dorewlts displeys the Ta- 
péacion operator in epherical coordinates. We perform e senarnute derive- 
tion here in order to Gsplay the transformation in a forum independent of 
representation. 


. 
"Th our notation 3 the upper index denotes the entire class of 
4 
Vesiables. Tor exemple, «” inolades Nea Ya» &xl %,, Gs well as 
& = 
r ge There ic no symbol y” 
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APPEADIN 2B 
SOLIPION OF "EIS CORRELATION AMPLITUDE BCUATTON 
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Figure 12. Coordinates for Long ‘ange Averaging Integral 
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We iliueteate the method used in constructing Fig. 6 by eveluat~- 
ing the pointe at « #@ 3.6 A. 


Ve first evnluete the cum Of 0 ee using the series expension 
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form of IV-{ 37) plotten in Pig. 5. 
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We consider the tera in IVik3) wizich gives the energy of the 
Plavinage Vertiusbatloa 
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m (3-4), the product over t is non~-@lagcnal so that tae may operate. 
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to lowest order in the energy, in G Beevais Inttice. Tet us consiaer 
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The fixet relationship stems fraa the fact that where @f1,2) has an 
eppreciable velue, near %5 * 0, g(1,2) is sowevhet more gently varying 
than eee, the gsecond is obtained oy inspecting Fig. 2. 
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there are four such pairs per atom, so that Bay ms 3,2° K, 

{no oatimating procediwe weed here demonstrates omly that H' is 
not the douinert term in the energy. The assumption thet H’ is omell 
end has little effect on the wave function must remain an esmeption of 
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In thie appendix ve obtain the factor 1/3 introdhoed into the 
energy in Eq. I¥-(27) by 0 wore detailed study. In Bq. TV-(27) the 
mateis elenont a 48, : (f, ft, ) = !t.) measures a “kinetic energy" 
of particle 1, using J .. > the pert of the veve function of particle 1 
which is determined ty corrcintions with perticle 3. the cum of all 
such contxlbutions is the cum total of the single particle energy recult- 


ing from the confining correlations. In methemetical symbols, 25 vf, 
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But in ow single-particle representation, sl1 the f's for 
ntoms eat the came distance a exo cphericaily syumetsic about the orl- 
gin, lattice site 1. The mibocsipte on the derivatives exe then super- 


fiuows. We can, therefore, write, for (G~2), 
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We can understand the factor 1/3 in enother woy, by applying 
CLaseicnl mechaxics and then identifying; the dynamical veslebles with 
the corresponding quantunme-nechanical operators. 
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We vioh to meagwre the kinetic energy of the peabicie 1 by ob- 


verving the lengths of 022 the u, ond teking a sultebdle average. ‘et 
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we agoguae thet we connot measure cur abschite position in space, but 
mevely ce measure the scalar quantities u “< (The problem ve wish to 
otudy is the inverse one, where wo measwume quantities in tems of 4, 
and deduce the corresponéing quantitien in use) 


We proceed by taking the tame derivetive of (G-5) and squaring: 
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Te kinetic energy is three times the avermge of the “kinetic energies” 
relative to the ?4xod points. 

in ou cnse, we esk for the apparent kinetic energy of relative 
motion to @ set of points j, in terms of the kinetic energy of a particle 


moving in m® poeth described by q, 
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By the correspondence principle, the right side of (G-G) represents the 
“effective kinetic energy’ of Eq. IV-(26), and the left side represents 
the modified kinetic energy of IvV-(27). We are, therefore, led to the 


factor 1/3. 
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fPPOTDIX 
WOUMALTOATION Of Til WAVE FUROLTON 


In thie appendin we attend to the detailed problem of normalising 
the grounl stete vave Sunetion 
in Chapter IV, we fond that the expression 


v(i) = -- 19, (291° 7 y(n) % tt 2%, 
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form of the },. 
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the trontment cf the sheli model, we cali upon thom to desorive the 
states pesmitted by the Pauli principle. Lach of these ptatves we will 
eomme to be ccoupied by tuo, one, or no particles. (Ve Imbel the states 
py £ oa got two particles per stete in our modified scheme, rather 
then by 3 * 2 : 1/2, for one particle per steta in the usual scheme.) 
Tren, since the particles ate kept in the vicinity of lattice sites, 

each weve muct be of a wavelength outfictent to apen two letiice sites. 


In Chapter IV, we wrote e general set of waves Tor m simple cubic 
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wouere the 6 4, wore chosen for each wave to put no nodes at inttice points, 
ex for n or 1/n off, to place meviua at lattice points, or for nor 
i/n even, maccima midway between lattice points. 

Ye will pow evaluate A in a particular case, where the lattice 
point %, is at the center of a soltd (RL a) by CH e) by (i, a). 
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We will estimate a for the case leading to the waves (1-3). 
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We cennot evaluate the swms in (R-G) wilese we know more about 
the reietive values of the Cc. me" we now show that, if we use the aseump- 
tionn preceding (1-2), (+6) tees on a particularly simple form because 


the C, 4k mre ail of the same magnitude. 
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Ve write the No body density matrix in a slightiy difverent 
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In (2-7), wo egoin neglect overlap. We, therefore, have 
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Because cach atete is occupied to the sme extent as any other stete, 
in the wey wo welte (1-7) an@ (1-3), all the pemautations give the some 
value. In (+o) we have written a narticular one, multiplied by the 
nudber of perminoations. 

We Jearned in Chapter IV that the correlations functions %, 
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we evaluated in Chapter iV as f,. 





ekp - 2 + tT HUN Rp 
“im. 1 + jet Lhe 
rT Se -  leeg +* 


(er me pe name ene ae 2 Pee of ee Oy me, 
om oc ong ee ee ee oe 
— tee er Le tien oe a GI ® | 
ae 
ee ee ee ee 
CUS" een ow ot we ere ie owe os «om 


nents let — en en peer lye - ~- 


° «7 —_—e UC seer @& 
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We, therefore, ean write (45) in the following form 
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Ve hove taxen the leeding texan in the Pa (The expansion we use 
in (H-~10) is expected to be valid, since oc 0° ¢ Will provide a cut off 
in the region of poor convergence. ) 

Ia the mum of (2-10), thore are atyye/3 temas with mn, ™ des 
ecm! 3 temas with ha 3, etc. This cawes about becouse the quentities 
Ya, > ifn 4? ifn, abe proportionnl to the eamponents of a vector in the 
Tesciproual lattice of our assimed eilaple cubic lattice. There are 
uy x u be cay? 3 Lattice pointe in the plane n, * (const) and an 
@gml nmugoer tn the rovtiected plone. 
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where wo have set the integral in the Gencetinater of (i-12) equal to A,» 
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As we did in (4), we can take all the cosines to be wiity. 
We then get: 
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Cdingle Particle Croétal 
. Beroarkes ond H. Primnkot?, Phys. Uev. 129, 960 (1960) 


The position acordinates of perticie 2, not inciwling 
opin (, inciudea opin) 


The coordinates of particle i, including spin by (2, 
does not inclule spin) 

Deexwees lelivin 

Angstrom unit 


Forms a Cleter detemzinent of vhat follows. The pesmutat: 
operator Yor Ferri perticles 


Yorms vo oyamotric pormriation of vhat foliows. 

Pesrutes ouly the coomtinuate infices 2, in the anti- 
oysmmetric gengc. 

Dirac § function 

The ‘wave number, of the single particle representation 


at) Q. 


vii) = conte 86, «= «g, © Gietance from Inttice site 1. 


fry) seg ar , 


r PRR: ‘as mx, i. wre x.) Density matrix of order (p) (See Chapter 


v(X'|K,) 


soe, 


B, 


Iv.) 
Density untrix of omer 2. (See Chapter IV.) 


Sie Cigpleacement between perticle i and perticle j. “his 
variable hes ea derivative 4/dx, but no derivative a/r, ;" 
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“ue anterpaerticle parmaoter, which is on independent vars- 
gbple in the Pluvinage formulation of tee Schrodinger 
equation. Jt has the constraint 
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Wich must be lmposed vhea evaluating matrix elementa can- 
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wo derivative O/fz. . 
we position of the center-of-umss of particles 4 and j. 
Gere is a derivative 3/Ox, wut no derivetive o/y, ; 
the correlation function, defined by formula ITit-(17). Tt 

. ral , 4 ry : oa = ae 
hao 2 series representation Flo/r, 4) O(r, j vy) » where 


@ is the step function, and a Yukava-type venvecentation, 
-y(r - My) 

The Qistance trom the Meld point k to the Inttice point 

The singlo-pecticle wave function, a function of speece and 

BDA0. 

She single-particle wave function, a Dertion of space only. 

i -— » 2 

che exror integral g{x) = (ofaleyy eat. 

7 £ tugs) where t is the position of an atom noar 
Stioe point +t. 


A series emprension for X, the correimtion function. 
An avbitrary function of (r, ry , 
ody 


fa, | - IR, 4| (See Pig. 4.) 
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The Gistance from the lmttice point denoting the single- 
pastiche xeference point to the shell of atoms at distance 


iL» @ tmetion of lattice geanetry. (See Mg. 4.) 
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